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ABSTRACT 


The present work deals with the steady state analysis and 
control of torque harmonics of the induction motor fed by three 
phase current source inverter. 

The analytical solutions for the rotor currents, stator 
voltages and torque waveform in the stationary dq frame 
attached to theSifrator have been obtained. The rotor currents 
and torque harmonic spectrums have also been obtained through 
frequency domain analysis. 

The next part of the thesis deals with the control of 
the torque harmonics by modulation of dc input current. A 
detailed study has been made of the modulation by exponential 
and cosinusoidal waveforms and their effects on various torque 
harmonics. A general nonlinear second order equation has been 
derived for obtaining the waveshape of the inverter input 
current required to produce a desired torque waveform. A 
general analytical solution of this equation has been done for 
the case of stationary rotor. Further for the case of rotating 
rotor and constant torque, this equation has been solved using 
phase plane analysis and numerical integration. The current 
profile for producing constant torque at any rotor speed in the 
motoring mode has been obtained. 



CHAPTER 1 


INTRODUCTION 

It is well known that an induction motor supplied by a 
current source inverter draws line currents with an approxi- 
mately quasi square waveform [l] . Consequently, the general 
analysis of the motor with sinusoidal voltages and currents is 
not directly applicable in this case. The most popular 
approaches for the analysis of ac motor drives, which employ 
non sinusoidal sources are harmonic superposition techniques 
using either symmetrical components [2] or method of multiple 
reference frames [3]. The analysis of the performance of the 
induction machine fed by a current source inverter has been 
carried out in [l] » [4] , [5] , [6] , [7] using dq axes model of 
induction motor. Several methods of solving dq equations 
have been proposed. In [4], [5], [6], the analysis of induction 
motor is carried out assuming ideal waveshape for stator current 
waveform. In [7] the analysis has been done using digital simu- 
lation. None of these give direct analytical expressions for 
calculating referred rotor currents stator voltages and torque 
waveform. In this thesis analytical expressions are obtained 
using d-q model of the machine for referred rotor currents, 
stator voltages and electromagnetic torque. Analysis has been 
done for both the cases of ideal current source inverter and 
the inverter with nonzero commutation time. For the latter, 
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the rise and fall of the stator currant during 
been taken of cosinusoidal nature [l]. 


commutation has 


It is possible to control the torque harmonics of an 
induction motor fed by a current source inverter either by 
modulating the dc current input of the inverter or by modula- 
tion within the inverter. In [8], [9] the effects on the 
torque harmonics by modulating the inverter input current by a 
particular exponential waveform has been studied. Modulation 
within the inverter is referred to as pulse width modulation ir. 
literature [lO]. Reference [10] gives the effects on torque 
harmonics because of the pulse width modulated current source 
invc3:ters. In this thesis a detailed study has been made of 
che modulation of inverter input current with exponential and 
cosinusoidal waveforms. It has been shown that the parameters t 
of the modulating waveforms control the torque harmonics. 
Optimum parameters for eliminating the dominant sixth harmonic 
have been calculated. A general nonlinear second order equa- 
xion has also been derived for obtaining the waveshape of the 
invertor input current required to produce a desired torque 
waver orm. A general analytical solution of this equation has 
been done for the case of stationary rotor. Further for the 
case of rotating rotor 'and constant torque this equation has 
been solved using phase piano analysis and numerical integra- 


tion. 



,The c’q model of the induction motor, with axes attached 
to the s to tor has been used throughout for the analysis. A 
brief description of dq axes model for three phase induction 
-no tor is given in Sec. 1.1. 

In Chaptc-r_2, the steady state analysis of three phase 
induction motor fed by an ideal three phase current source 
inverter has been done to compute the referred rotor current 
and the stator voltage waveforms. Tho analytical solutions are 
obtained for rotor currents and stator voltage for three inter- 
vals of the inverter period. Boundary conditions are verified 
analytically for the c'use of stationary rotor. The resulting 
expressions for the case of the rotating rotor are complicated 
and hence, closed form equations have not been obtained. A 
computer program has boon written to compute tho rotor currents, 
/analytical solution results in a considerable reduction of com- 
puter time compared to the numerical solution of dq equations. 
The boundary conditions for the case of rotating rotor arc 
verified with the computer program. The values of the rotor 
current harmonics have also been obtained through frequency 
domain analysis of the induction motor. It has been shown that 
the rotor harmonics as obtained from tho analytical solution 
ore in close agreement to those obtained through frequency 
domain analysis . 

In Chapter 3, the analysis similar to that in Chapter 2 
has been done. In this case tho invertor is assumed to have 
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a nonzero commutation time. During this commutation time., 
the rise and fall of the current is approximated by a 
consinusoidal function [ 1 ] . 

In Chapter 4_„ the electromagnetic torque of the current 
fed induction motor has been computed. Two methods to compute 
the torque have been given. One obtains the analytical solu- 
tion of torque using the time domain expressions of stator and 
referred rotor currents in dq frame. The other computes the 
torque harmonic spectrum using the harmonic components of the 
stator and rotor currents in d*-q frame. Using these procedures - 
the torque- waveform and spectrum has been studied for tho case 

of ideal current source invertor supplying an induction motor, 
to the -adulation of the inverter input dc current by exponen- 
tial noeulstion. 

in Che ’ter_5, a detailed study has been made of tho 
effects on torque harmonics due to tho modulation of the inver- 
tor input dc current by exponential and cosinusoidal waveforms. 
Tho exponential modulation for stationary rotor case gives con* 
stant torque. Cosinusoidal type of waveform may be present in 
the invertor input current due to imperfect filtering in the 
dc link. It has been shown that it is not possible to get an 
exactly constant torque with these* typos of modulation for^the 
case of rotating rotor. However, by a certain choice of the 
parameters of the modulating waveforms, it is possible to 
reduce certain torque harmonic. 
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Chapter 6 deals with the calculation of stator current 
profile of the induction motor fed by current source inverter 
to produce a specified torque waveform. A general nonlinear 
second order equation has been derived for obtaining the 
v/aveshape of the inverter input current required to produce 
the desired torque waveform. A general analytical solution of 
this equation has boon done for the caso of stationary rotor. 
Further for the caso of rotating rotor and constant torque 
this equation has been solved using phase plane analysis and 
numerical integration. The current profile for obtaining a 
constant torque at any rotor speed in the motoring mode has 
been obtained. 

1.1 DQ MODEL FOR INDUCTION MOTOR 

The development of the dq model for the induction motor 
using arbitrary revolving frame has boon done by Krause [ 3 j . 
In this section the equations of motor performance and the 
transformations from three phase balanced system to the dq 
frame have boon outlined briefly using [3]. 

Fig. 1.1 shows the angular relation of the stator and 
rotor axis of a three phase machine with the third set which 
is an orthogonal set (dq axis) rotating at an arbitary ele- 
ctrical angular speed d©/dt. It is clear that a^-b^-c-^ 
is fixed in the stator and a2“b 0 -c 9 sot is fixed in rotor 



1.6 


and hence rotates at an angular velocity to . Subscript 1 is 
used for stator quantities and subscript 2 for rotor quanti- 
ties. The time zero angular relationship between three sot 
of axes can be- selected arbitarily? however, it is convenient 
to assume that at time zero ,a^ end ’d s axes coincide. 

The transformation equations which for the balanced systc 
can be correlated to the angular relation of the axes shown in 


Fig. 1.1 ; 

1TQ 

written 

as follows 

O 

O 



f dl = 

2/3 

[f 

al cos0+ 

f bl ' 

cos(0 - 

■§*>♦*01 c 

os(© + ) ] 

( 1 . 1 ) 

f , = 
ql 

2/3 

[~ 

f al sine “ 

f bl 

sin(0 

2ll\ r- 

~ T”' cl 

sin(©+ ^L) ] 

(1.2) 

-p _ 

L d2 ~ 

2/3 

[f 

a2 cos|3 + 

f b2 

cos( p~ 

‘ §“> + f c2 

cos( j3+ I 2 -)] 

(1.3) 

V = 

2/3 

[- 

f a2 Sinp- 

f b2 

sin( p~ 

■ I 2 ') “ f c2 

sin(j3+ §~)] 

(1.4) 

where 














P = 

0-0 
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(1.5) 


In those equations the variable f can represent either voltage, 
current or flux linkage. The equations arc restricted in that 
the instantaneous angular displacement 0 of the arbitary 
reference frame must bo continuous finite function. 

It has been shown in [ 4 J that the equations of motor 
performance for squirrel cage induction motor for a stationary 
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reference frame are 
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(1.7) 


( 1 . 8 ) 


where m is the number of phases, P is the number of poles, L-^ 
is the leakage inductance of the rotor, 'the leakage 

inductance of the rotor referred to the stator, M is mutual 
inductance and ' p ! refers to 'd/dt'. In these * T * is the 
instantaneous torque of the motor. 

In order to solve the above dq equations, the three phase 
input excitation currents are transformed Into two phase dq 
currents. The transformation connecting 3 phase to 2 phase 
variables with assumed directioh of coordinates as shown in 
Fig. 1.1, for tho case of balanced stator current gives 


1 dl 1 al 

1 ql y 3 ' x bl cl y 


(1.9) 


( 1 . 10 ) 
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Ficj. 1.1 


Axes of throe- phaso symnotrical machine 
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CHAPTER 2 


STEADY STATE ANALYSIS OF THREE PHASE INDUCTION MOTOR 
FED BY AN IDEAL THREE PHASE CURRENT SOURCE 

INVERTER 


2.1 INTRODUCTION 

The three phase squirrel cage induction motor fed by an 
ideal square wave current has been analysed in this chapter. 

Fig. 2.1 shows the waveforms of the line current in the motor 
produced by ideal current source inverter. The dq model of 
the induction motor, with axis attached to the stator (Sec. 1.1), 
is used for the analysis. In the* analysis, the steady state 
solution has been obtained. This solution is for a constant 
invertor frequency and a constant rotor speed. The inverter 
frequency is adjustable by tho inverter gating circuits. 

The three phase current is transformed into the d-q 
current using tho equations (1.9) and (1.10). The transformed 
d-q currents, for all the six intervals of the inverter period 
are listed in Table 2.1. In Table 2.1, 

I = 2/3 I dc (2.1) 

whore I d is the dc link current fed to the inverter. 

Tho waveform of the d-q currents is given in the Fig. 2.1. 
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The waveform of the rotor current is obtained by the 
solution of equations (1,6). In the case of current fed 
induction motor, the above equations are decoupled from eqn. 
(1.7) which correspond to the phase voltage of the motor. This 
Is in contrast with the voltage fed induction motor, in which 
case all the four equations of the matrix equations (1.6) and 
(1.7) have to bo solved simultaneously. 

Since wo are solving for rotor currents at constant rotor 
speed, the problem can be viewed as obtaining the steady state 
solution of the differential equations with constant coeffi- 
cients. This has boon done using two well known techniques 2 

(1) Steady state solution in time domain determined by 
obtaining homogeneous and particular solutions and using two 
boundary value conditions. 

(2) Steady state solution in frequency domain obtained by 
considering the input as the summation of sinusoidal signals, 
through fouricr series of stator current and calculation of the 
transfer function for those inputs from the equation (1.6). 


Equation (1.7) has been used in this chapter to obtain 
the voltage waveform in the d-q frame at the stator terminals. 
These equations have terms containing 'P^dl* anci '^ql 1 * 

1 idi’ and ’ i are discontinuous (Fig. 2,1) the voltage wave- 
form is shown to have impulse. 
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Section 2,2 deals with the calculation of the rotor 
currents, using the time domain analysis. In this the eases 
of the stationary and the revolving rotors are dealt with 
separately. In case of the stationary rotor, the motor 
equations relating the dq s to tor and rotor currents, from 
equation (1.6) turn out to be first order and decoupled. The 
discontinuities in the rotor currents duo to discontinuities 
in the stator current can be computed from equation (1.6). 


In case of revolving rotor, equation (1.6) has two 
coupled equations for the rotor dq currents, i^p and ip. A 
single equation for i^p or iqp can k® obtained from equation 
(1.6) by elimination. However, the resulting single variable 
equation is a second order equation which contains the terms 


as 'p 2 i dl ' and 'p 2 !^ 1 . Since the stator. dq currents, i 
i ^ are discontinuous, these terms correspond to the deriva- 
tives of the impulses. To avoid such terms, the two now 
variables Xj and X 9 , have been defined instead of the rotor 
current variables, i d2 and iq 2 « Those variables and Xp 
have boon shown to bo continuous and arc proportional to the 
flux linkage in the machine. Those are sometimes referred to 
as pseudo rotor currents in the literature [4]. 


dl 


and 


The solution of the equation (1.6) in time domain 
requires two boundary conditions. These have been obtained by 
noting the fact that under steady state at each multiple of 
60 J in rnt, to being the invertor frequency^ is similar. The 



stator mmf jumps in space by 60° due to the particular form 
of i ,i,_ and i (Fig. 2.1) which have discontinuities . The 
stator mmf instead of rotating smoothly as for sine wave 
currents jumps by 60° in space for every T/6 period of the 
inverter , T being the inverter time period. Thus, under 
steady state the rotor mmf must also move by 60° in space in 
each T/6 interval so that each interval is similar. The 
resulting boundary value relations based on above arguments 
are given in equations (2.5) and (2.6), 

In Section 2 .2, the solution for the initial value of 
rotor currents has been obtained independently for each inter- 
val. From these it can bo seen the current solution obtained 
at the start of each interval is same as that obtained for the 
end of the previous interval. Thus it is sufficient to solve 
for the initial value of rotor current in any one of the six 
intervals . 

Section 2.3 deals with the technique of frequency domain 
solution. For the rotor currents in this scheme. the expre- 
ssions for the rotor harmonic currents have been obtained in 
terms of the harmonic components of the stator current. The 
cases of stationary and rotating rotor have been dealt with 
in one general expression. 

In Section 2.4 , the voltage expressions in dq frame at the 
stator terminals have been obtained in time domain. The results 
of Section 2.2 and matrix equation (1.7) have been used for the 
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c .lc ulcer on 


Foi 


verification of the solutions through 


ho tli the techni : ucs for rotor currents, two computer program 
hive boon developed. The first program computes the time 
dome. in solution of the rotor current and then obtains the 
harmonic components of this solution. The second programs 
directly computes the har ionic components of the rotor 
current by the frequency domain technique (Section 2.3). In 
Section _2.5 the case of a particular induction motor has been 
studied. Using the parameters of this induction motor, it has 
b-'en shown that results obtained through both programs ere in 


: lo: 


?roemcn-c, 


This verifies these ms u Its. 


2.2 ANALYSIS Or INOUCTIOM MOTOR FOR ROTOR CURRENTS THROUGH 
TL" ic DOMAIN 


In this section, the steady state solutions through the 
time domain technique, Section 2.1, for the rotor currents 
have been obtained. For this the motor performance equation 


(1.6) has boon solved. The c'ses of stationary and the 
rotating rotor have been dealt with separately. The equation 
(1.6) has been solved for the first three intervals (I, II, III) 
o < cot < it. The stator current waveform (Fig. 2.1) is 
symmetric about cot •= %. Therefore, solutions for intervals, 
IV, V and VI, it < cot < 2it,are negative of the solutions of the 
intervals I, II and III respectively. 
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For the analysis, each interval has been associated with 
the discontinuity at the start of the interval, Vfhile solving 
equation (1.6), each interval has been divided, into two sub- 
intervals. One corresponds to the point of discontinuity at th 
start of the interval and the other to the remaining period of 
that interval. 

The following symbols have been used for the analysis. 


bo = Wh = o") ; l q0 = = 0-) 


I d2 ” J 'd2' 'i 


Uo(t, = 0 + ) r X ' 2 = i q2 (t ± = oh 


U,(t, = TJ ; I£ = i q2 (t ± = T c ) 


<32 “ <32' i “ ‘c 


( 2 . 2 ) 


where T is the duration of each interval, i.e. (T = Tr/3(o) . 

V V 

The time ’t^’ is the time within the interval where i denotes 
the interval number i.e. i = 1,2,3,. Thus tg , ±2 and tg are 
defined for intervals I, II and III respectively and in these 
we have (t^ = t) , (tq = ( t—ix / 3to ) ) and (tg = ( t“2n/3(j ) ) . As 
already mentioned that each interval is associated with the 
discontinuity of stator current at the start of that interval. 
Thus ( t^ =0*) refers to time just before the discontinuity 
has been considered in interval ’ i’ and ( t^ - 0 + ) corresponds 
to the time just after the discontinuity has been accounted for. 


2.2.1 Case of Stationary Rotor 

From equation (1.6), the equations relating the dq stator 
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and rotor currents for the stationary rotor (10 

+ (V 2 + ^22p ) ^c!2 =: ^ 

ii.lpi.ql + ( Y '2 +L 22 p ^ ZL q2 = 0 

In this sub-section we are first obtaining the solution of i 2 
and i r , 9 from above equations with the arbitary initial values 
of I and Then, tho fact that steady state dq currents 

at tho e 'id of each invertor period make the airgap rnmf wave 
turn tc/ 3 electrical radians, is used to compute I^ 0 and I . 

The values of dq currents at the end of the inverter 
period are related to the initial values by the following 
equations [4], [ll] 

132 = U d0 - V3 Iq 0 )/- (2.5) 

^2 = dqo + 'f 3 W /2 < 2 - 6 ’ 

where I^d, and are as defined in equation (2.2). 

In Section 2. 2. 1.1, the rotor currents expressions for 
interval I have been computed. Sections 2. 2. 1.2 and 2.2. 1.3 
give the rotor current expressions for intervals II and III 
respectively, which can be obtained by analysis as done in 
2. 2. 1.1, Sec. 2. 2. 1.4 shows that it is sufficient to calculate 
I ( j 0 and I for one interval only. 


= 0) are, 

(2.3) 

(2.4) 
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2 . 2 . 1.1 Interval I (o < cot < tz/3) 

As t^ = t in this interval, the time is denoted as t 
i n th is interval. 

C aJL c u i_a i t i o _n of_ currents _at t_ =_ 0 

From Table 2.1, at t = 0, 

Ai dl (t=o) = i dl (t=0 + )“i dl (t=0“) =0 (2.7) 

Ai ql (t=o) = i cql (t=0 + )- i c , 1 (t= 0 “) = -*f3I (2.8) 

v ^ jl. 

Integrating equation (2.3) from t = 0 to t = O', 

0 + 0 + 0 + 

M J (pi dl )dt + r 2 / i d2 dt + L 22 J (pi d2 )dt = 0 

t=0" t=0" t=0“ 

Since i d2 is finite, this gives 

M[idi(t= 0+ )“i dl (t=0'“)]+L 22 [i d2 (t=0' f )- i d 2 (t= 0 “)] = 0 

So frora (2.7) 

J,lii d i (t=0) + L 22 Ai d2 (t=0) = 0 (2.9) 

So from equations (2.2), (2.7) and (2.9) 


^ 6.2 *do 


( 2 . 10 ) 
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Similarly from (2.4) using (2.4) oncl (2.3) 


i^ 9 (t=0) = 


iu 


-22 


•ql 


So, 



+ f3 Kj_ 


v/'nere 


1V 1 


l 1 ,-, 

L 09 


-alculntion of currents for 


From Table 2.1, during 


t >, 0 

this period 


i dl = -31/2 

iql = ■''/'3I/2 

Solving (2.3) and (2.4) using (2.14) and (2.15) 
i r - ; 2 ( t ) = 1^2 e"p(-at) 
i c 2 ( t ) = X q2 ex 'P(“at) 


where 


a 


r 9 /L 


22 


( 2 . 11 ) 


( 2 . 12 ) 


(2.13) 


(2.14) 

(2.15) 

(2.16) 
(2.17) 

(2.13) 


and I^o and I^ 9 


are as defined in equations (2.2). 
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Substituting (t = te/3oj) in equations (2.16) and (2.17) 
ad using (2.10) and (2.12) we have 


T 1 ? 

1 d2 


I do exp(-aT c ) 


(2.19) 


I M 0 = (1 +V3IO exp ( --al ) 
q.i ' qo * 1 1 ' c 


( 2 . 20 ) 


Solving for I, and 1^ from equations (2.5), (2.6), (2.19) 


md (2.20) 


qo 


'00 


i f ) i- ’ 

3 lC 1 "'2 

p' ‘ '“'Tn* 


( 2 . 21 ) 


qo 


-V3 k-j kA 

■ gb [exp(~aT c ) 


1 / 2 ) 


( 2 . 22 ) 


vi ne ro 


k 2 


exp('*aT ) 


(2.23) 


exp ( -2 aT ) -exp ( -aT_ ) +1 


(2.24) 


2. 2. 1.2 Interval II (i.e. | < art < ™ or o < t 2 < §) 


The rotor current results have been obtained for this 


interval preceding similar to Sec. 2. 2. 1.1, 



2 all 


At t 2 = 0 


' t — T 

'd.2 go 


P i c » 
2 1 


■\f O 

I'., = I + - J 7v k* 
c|2 qo 2 1 


For t 2 > 0 


(t 0 ) =1^ exp(-at 0 ) 


d2 v 2 • 


(2.25) 


i ( to ) = 1^2 


( 2 . 26 ) 


a no 


^r 3 


K 


I do 


O iV -t 

= | - ,1- q - 2 - [exp(-aTj] 


(2.27) 


I 


qo 


V"3 Ki XJ> r exp(-*aT ) 
~~ r — [1 - --y---] 


(2.28) 


2% 


2.1.3 Interval 111 (i.e. < t < % or o < tg < •§■) 


Similarly, for this interval. 


at t 0 = 0 


I d2 I do “ ’2 K 1 


X » = i «. )T3 k , 

1 q2 qo 21 


(2.29) 
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for t 0 > 0 

O 


i d2 (t 3 ) = X d2 ex P(" at 3 ) 


(2.30) 


1 q2 (t 3 ) 


I q2 ex P(“ at 3 ) 


(2.31) 


an: 


'do 


qo 


K* ’ 

3 1 2 


C, 


[exp(-aT )~l)] 


*.r. 


• K o 
g - c * Cexp(-aT c )+l] 


(2.32) 


2. 2. 1.4 Discussion 

The rotor current expressions for intervals IV, V and VI 

vnitefvo\s 

can be observed directly from those in S ee-t i rQ - n-c I, II and III 
respectively because of symmetry property. He also observe 
that for this case of stationary rotor, the closed form 
expressions for L. and I _ are obtained. In the above 

C i Q Cj O 

sections, the analytical expressions of I^ Q and I have been 
computed independently for each interval. The expressions 
of the final values, at the end of each interval, can be 
obtained from equations (2.3) and (2.6) using the expressions 


of I 


do 


and I 


qo 


obtained for that interval. These values are 


summarised in Table 2.2 


It is seen from Tables 2.2 that the rotor currents at 
the start of each interval, i.e. at (t^ =0) is same as that 
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obtained at the earl of previous interval i.e. at = T c ). 

Thus ’vc need to solve for 1^. and I _ for one interval only. 

GO Cf O 

2.2.2 Case of Rotating Rotor 

From equation (l.o), the equations relating the dq stator 
e nr rotor cormts with co, as the angular velocity of the 
■ /cor in electrical radians per second are 


:uo 


r x dl 


'M pi^ t ^ "p^^ppP ) i r /9 '* ^-’op^/pi 


2 22 1 ' e2 22 r c!2 


(2.33) 


!'u 


>x d 3 + kco r ~c, 1 ' L 22"°r 1 c:2 + ^ r 2 +L 22 p ^d2 “ 


0 


(2.34) 


These equ' tions ore coupled in contrast with the stationary 
ro case, equations (2.3) and (2.4) where these are decoupled. 
Further pi cll and pi^-, are impulse functions as i^ and i ^ ar s 
discontinuous. The solution due to impulse functions has been 
■: v. -0 fallows ", 


To obtain the changes in rotor currents at any time instant 

due to change in stator current we integrate equation 

(2.33) from t = t rt to t = t* . This gives 

o o 


o 


4 . 

t; 

•4 




v) r i dl dt + M f Pi q l dt + 


t=t 


o 


t 

o 

r» 

r 2 J 
t=t" 


t 


i q2 dt + l 22 


f p i,^dt 


t=t 






L 0n J tO ^ \ d '>«. 0 


(2.37) 


o 
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ir rents are finite, 


It = 0 wht 


q 2? 


o.-, or cl. 


(2.30) 


■s , from (2.37), (2.30) 


.[i .(t=t + ) 

u Csl o 


V u=t 4 ]+ 42 [ 4?. (t= v 


J- o V O— O 

a 2 o 


=t")3 = o 


(2.39) 


1 o2( t =t 0 ) = 


f ' \ 

■ql v V 


(2.40) 


where ( t-t ) corresponds to any arbitary time instant and ’ 
corresponds to a jump in the current variable. Similarly, 
from equation (2.34) 


<A i d2^ t “’ t o^ ~ 




(2.41) 


Thus the effect of the discontinuity in i^ and i ^ at 
tj, = 0 is a corresponding discontinuity in iq 2 and i 2 , given 
by equations (2.41) and (2.40) respectively. This is equivalent 
to stating that the response due to 'pi^j_ ! and * Piqi' i m P u l ses 
in equations (2.33) and (2.34) can be looked upon as sudden 


changes in the initial conditions of i^ 2 and i 2 < 


Therefore, 


.t is possible to solve (2.33) anct (2.34) as a set of simul- 


For the interval t^ > o, because i^ and 


tanoous equations. 
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"ql 
a s 


are constants (Fig. 2.1) and pi^, P^qj_ are zeros. Thus 
olution similar to Sec. 2.2 of stationary rotor can be 


9 


obtained . 


However, equations (2.33) and (2.34) can be solved more 
conveniently if a change of variable is done from i^ 2 and i^ 2 
to X 9 and X 1 respectively. These variables, >0, and are 
defined as, 


X 1 ~ 1 q2 + L 22 


(2.42) 


x 2 = l d2 + 


.JL i 

L 2 2 d 1 


(2.43) 


It is shown in the following that these new 
variables and X 2 are continuous functions of time. 

The change in at a arbitary point (t = t ) can be 
obtained from equation (2.42) as 


£ X 1 (t=t ri ) =^i q2 (t=t Q )+ i^i„i(t=tj 


J 22 


ql' 


(2.44) 


Using (2.40) in (2.44) we get 


A Xj_( t=t Q ) = 0 (2.45) 

Similarly from (2.43) and (2.41) 

A.X 2 (t=t Q ) = 0 (2.46) 

These new variables X^ and X 2 are referred to as pseudo- rotor 
currents [4] and are proportional to the rotor flux in the 
machine . 
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In equations (2.33) and (2.34), eliminating Soar ’ pi^’ 
and , piq 2 * by using (2.42) and (2.43), we obtain 


P X 1 = w r X 2 " ai q2 


pX 2 = -u r X 1 - ai d2 


where 


- 2 /ij 22 


(2.47) 

(2.48) 

(2.49) 


Substituting for i^ 2 and i^ 9 in equation (2.48) and (2.49) from 
(2.4^) and (2.4^) gives 


L 1 “ r*' v 2 


aX i + r 


(2.50) 


pX 2 =-(o X, 


aX 2 + L 22 ai dl 


(2.51) 


5 © 7 

Dif forontiating (2.4*7) and substituting for pX 0 from (2.44) and 


similarly differentiating (2.51) and substituting for pX^ from 
equation (2.48), gives 


2, . 2 W 


where 


= Kl^ql+Wdl+KlPiql 

(2.52) 

= K 1 ai dl -K 1 io r i ql +K 1 pi cll 

(2.53) 

— - Ma/L 22 

(2.53a) 

and (2.53) are the final equations 

rotor current expressions. 

which are 
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2. 2. 2.1 Interval - I ( o < cot < %/3 ) 

.’.t tot = 0, from continuous property of X-, 9 %2 

2 = 0 A X ? _ = 0 (2.54) 

For o < tot < it/3, from Table 2.1 


V 

il 

1 — 1 
*D 

*rH 

-31/2 


± ql ” ' 

'731/2 

(2.55) 

From (2,52), for 

this interval 



9 9 9 K I a/" 3 1 K. -i to 3 r 

(p +2ap+a +oo^) X-j = — ' ~ " 2 ”“ ' ** ■*-— vs--' — 


The solution of above is 


Xj(t) = exp(-at) sin oo^t + cos to^t ) ] - (2.56) 


V'/herc 


^4 


y3 Sw^+a) 


'2 ‘7 ' 2 ’ "“ 2 \ 
(a +toJ 


(2.57) 


ant; 


C 7 ., D; are constants, 
oi ? ol 


/: io 

•,r 3 

= ^i! 

x i 

! 

t i =0 

J 

V 

" v 20 

— %/■ 

- -' v 2 

; 

t - =0 

1 1 


• 

' C 10 

cl X . 
= dt~ 

L 

‘ V 0 

« 

J 

• 

,( 20 

dX 2 ( 

“ eft ' 1 1 . = c 

(2.58) 

X” 

/V 1 

_ - a j 

- ! 
-L ; 

t,=T 

• 

9 

Y «* 

2 

— x 0 

t,= T„ 
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where i corresponds to interval. Here i = 1 and t^ = t. 
Using above in (2.56) gives 


D ol = X 10 + K 4 

(2.59) 

C 'ol = (X 10 +aD il )/w r 

(2.60) 


Solving (2.50) at t = o, using (2.55) and (2.58) 

*10 = to r X 20 “ aX 10 " (2.61) 

Similarly solving (2,53) for Xq gives 

X 2 (t) = exp(-at) [C^ sin co^t+D^ cos w r t] + *6 (2.62) 


whe re 


= 


V3 

2 


I^Kuy- 

(a 2 +a)p 


a’f3) 


D o2 “ *20 “ K 6 
C o2 = ^*20 +aD o2^/ w r 


X, 


20 


-co X , 


r 10 


“ aX 20 " 'S K 1 I 


Computing X^ at t = T c from (2.62), using (2.58) 


J\ i 


t=T 


X 1 = C il Y 3 + D ol y 4 


K 


4 


where 


Y, = exp(-aT ) sin co T 
O C X* c 

Y 4 = oxp(~jjT c ) cos (0 r T c 


(2.63) 

(2.64) 

(2.65) 

( 2 . 66 ) 

(2.67) 

( 2 . 68 ) 


(2.69) 
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Substituting for 0^, in (2.67) from equations (2.59) 

and (2.6C). Using (2.61), we have 


wnc ro 


where 


X 1 = Y 4 X 10 + Y 3 X 20 + Z 3 


Z 3 = C 5 Y 3 " W4 


= [- Kl I + K 4 a]/w 3 


(2.70) 

(2.71) 

(2.72) 


Similarly computing at t = T and solving as above we get 

(2.73) 


y n — „y y .’ r v Y -+- 7 ’ 
2 ~ y 3 10 Y 4 20 4 


where 


where 


Z 


4 


r ? 
u 6 


C 6 Y 3 ” K 6 Y 4 + K 6 


[- 1 K X I - aK*]/o>, 


"6 J/LU r 

writing equations (2.70) and (2.73) in matrix form 


(2.74) 


(2.75) 


i s 


— 



— ~i 


— 

! v ? ? 1 

^2 ' | 

j 

1 

! 

i 

- Y 3 


■v 

*20 

+ 

7? 

i X 1 ! 1 
! l j 


Y 3 

Y 4 


. ""10 


7 1 

j » 


| 

- 


J 


— — 


(2.76) 


From the boundary condition property, the relation between 
the initial and final values of current of an' invertor 


interval, 

I X 2 


x r 


1/2 

V3/2 


~V"3/2 

1/2 




X, 


20 


v 

A 


10 


(2.77) 
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From equations (2.76) and (2.11), solving for 


X 20 


- Y 4 ) 

-( f l - v 3 )" 

-1 

N 

4*- 

) 

X 10 


<% y 3 ) 

- Y4) 


7 » 


X 20 and X 1Q are known from (2.78) after calculation of constan 
in the order K py K 4 , K 6 , Y 3 , Y 4 , C”, C^, 2^ and Z^. 

The pseudo currents (t) and X 2 (t) are now known from 
equations (2.56) and (2.62). The actual rotor currents i ^ 2 an ~ 
iq2 can be obtained from equations (2.44) and (2.45)* 

2. 2. 2. 2 Interval II (i.e. ■§ < cot < §2 or o < t 2 < g) 

Solving as interval I we obtain 


X x (t 2 ) = exp(~at 2 ) [cos 

sin ‘"H cos 

(2.79) 

where 



T3 K, la 

7 

(a +m r ) 


(2.80) 

D o3 = X 10 + ^ 


(2.81) 

C o3 = X 10 + aD o3 


( 2 . 82 ) 

and 




X 2 (t 2 ) = exp(*-at 2 ) [Co^ sin w r ' t 2 + ^o4 cos w r ^2^ + ^10 


( 2 . 83 ) 


where 

K 


10 


V3I 

(a +w r ) 


( 2 . 84 ) 
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D 04 X 20 " K 10 


^04 ~ + a ^n./i]A° 


20 


04- 


(2.35) 

( 2 . 86 ) 


10 » ^0 5 A io cUiU A 20 corres P onCt x0 


It should be noted that X 1Q , X-,^, X, ^ and v 
definition of equation (2.58) with i = 2. 

For initial condition of this interval we have the 
relation 


'r 1 

' v 20 

! 


(-2- - V 4 ) -Of - Y 3 ) 

-1 

r z> "I 
8 | 

i 

o 

1 — 1 

>< 

: * i 

. ( ir | - Y 3 ) (| - Y 4 ) 


! z- 

L 7 J 


(2.87) 


where 


z 6 

== 

C i2 Y 3 “ K 10 Y 4 

+ K 10 


(2.88) 

Z 7 

= 

C il Y 3 + K 7 Y 4 

" K 7 


(2.89) 

C l2 

— 

“ K 10 a/w r 



(2.90) 

C il 

= 

[-*V*3 K X I + K ? a] 

Ao r 


(2.91) 

X 

o 

and 

X 2 Q are solved, 

the pseudo 

currents Xj_(t) 

and 


Xq(t) are defining over the entire interval tc/ 3 to 2n/3. The 
actual rotor currents can be found out using equations (2.44) 
and (2.45) . 
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M-p m* 

2. 2. 2. 3 Interval III (i.e. 3 ” 5, < 71 or ^ < ■j) 

Solving as for interval I we obtain for this interval. 
X 1 (t) = exp( -atg) [cos ' sin u r t^ + cos 0 J r t 3 ] + K ^ 3 (2.92 


x 2 (t) 

= exp(»at 3 ) [Co£ sin w r t 3 + D ^ 6 cos o> r t 3 ] + K 16 

(2.93 

where 



Jc 3 

= t - 2%/3 

(2.94) 

K 13 

f3K x I (V3w r -a) 

_ — 2 — - 

( a +w r ) 

(2.95) 

K 16 

V3K.I (w 4f3a) 

_ — 5 — 

(2.96) 

D o5 

= X 10 “ K 13 

(2.97) 

D o5 

= [X 10 + aD o 5 ]/w r 

(2.98) 

D o6 

= ^0 " K 16 

(2.99) 

C o 6 

= [X, 0 + aD o6 ]/“ r 

(2.10C ' 


For the initial conditions, the following relation is obtained 


X 20 


(? ~ Y 4^ 

~( r i - y 3>" 

“1 

~ Z i 2 ~ 

x io • 

j-N. 

1 

Si - v 3 > 

<2 - Y 4> 


Z Il 


( 2 . 101 ) 
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where 


Z l2 


33 C i8 Y 3 - K 16 Y 4 + K 16 


z il ~ C 17 Y 3 “ K 13 Y 4 + K 13 


where 


r V3K1I 


*13 J r 


Ol/ T 

C ls = [ - K 16 a]/« a 


( 2 . 102 ) 

(2.103) 

(2.104) 

(2.105) 


Thus, the actual rotor currents have been found out for this 
interval. In rotating rotor case it can be seen that analytical 
expressions are complicated and closed form expressions are not 
obtained for dq rotor currents. A computer program has been 
developed to calculate initial and final values of currents 
for these intervals in rotating rotor case. The boundary condi- 
tions have been shown- to be matching for the case of induction 
motor in Section 2.5. 


2.3 ANALYSIS OF INDUCTION MOTOR FOR ROTOR CURRENTS THROUGH 
FREQUENCY DOMAIN 

The expressions for different harmonics of the rotor 
currents have been computed in this section. The cases of 
stationary and rotating rotors have been dealt with in one 
general expression. 
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The analysis proceeds in two s teps.F irs tly , is the 
calculation of the harmonic components of the stator current. 
Then equations (2.52) and (2.53) are used to compute the 
harmonics of pseudo rotor currents and With this 

result, using equations (2.42) and (2.43), the harmonics of 
rotor currents can be obtained. 

For calculation of the harmonic components of stator 
current, i^ and i ^ (Fig. 2.2) are resolved into Fourier 
series as 

oo 

i ,, = E a cos nut (2.106) 

ox n=l n d 


i , = E b sin nut 
qx n=l n q 


where 


for all even values of n 


a =0 
n d 


b n =0 
ci 


and for all odd values of n 


61 . nit 

1 = -.i sin iAL. 

n^ nit 3 


b . . 4f|X sln 2 ™ 

n nit 3 

q 


(2.107) 


(2.108) 

(2.109) 


( 2 . 110 ) 


( 2 . 111 ) 
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It can be seen from (2.110) and (2.111) that 
triplensi.e. n equal to a multiple of 3, a 

n d 


Thus 


for all 

and b = 0. 
n 


f or n 


a 


n 


d 


b 


n 


q 


1,5,7,11,13,17, ... 

(2.112) 

- -nV '"I 

rax 3 

(2.113) 

_ 3£l t 
rax 

(2.114) 


and all other harmonics are absent. 


In the remainder of this section only o dd and non triplent 
h armonic s are considered from (2.106) and (2.107) 


. 61 . rax t 

1 d 1 = 2 - "n% sin cos nwt 

1 n 


i = E - sin nut 

Q n nTc 

*1 n 


Writing for mth harmonic of i , and i 

c\ 1 q 1 


bl = cos (™t + 0! md ) 

m 


1 ql m = I mq cos ( m “ t + % q ) 


So, from (2.115) and (2.116) 


(2.115) 


(2.116) 


(2.117) 

(2.118) 
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w 

= l X mq» = ii m i 

(2.119) 

a md 

=5 0 

(2.120) 

a 

rnq 

= - %/2 for m = 1,7,13,19, ... 

(2.121) 

a 

rnq 

= + tc/ 2 for m = 5,11,17,23,... 

(2.122) 


We define a variable * p' as 


1 if m = 1,7,13,19, ... 

P = (2.123) 

0 if m - 5,11,17,23,... 

then from (2.117) and (2.118) 


i dlm = cos rnwt 

i qlm = I m cos + ^/ 2 ) 


(2.124) 

(2.125) 


wh e re 

^ - ("D P W ■ a nd < 2 - 126 ) 

Solving (2.52) for rath harmonic of X^ denoted by X^ m , equation 
becomes 


( p^+Pap-i-a^+co 2 ) X lm 


K. ai , + K, to i + K, pi . 

1 qlm 1 r dim 1 H qlm 


This equation can be written in phasor form as 



where 


A t 

r / 2 2 2 2v2 <.2-1 1/2 

= L(~m w +a +u> ) +(2aa)m) J 

(2.128) 


= tan” 1 [2anxa/(a 2 +w^-m 2 w 2 )] 

(2.129) 

■^qlm 

= ^ Z.bD p 1 

f2 z 

(2.130) 

^dlm 

= is /o° 

V2 ■' ' 

(2.131) 


and X lm corresponds to the phasor of X lrn (t) from (2.127), 
(2.130) and (2.131) 

K,a -o, K,u 

X lm - <A7 /=*r) X qlm + t~xf ^ X dl m + 
mwK, 

+ (-^-~ /j 90 - a T ) I qlm (2.132) 

— ,> 

This can be visualised as if X,„„ is obtained from summation 

lm 

of three phasors which can be obtained from the I q ^ m and I^m 
by passing these through appropTi ate amplifier, as shown in 
Fig. 2.2. 

To obtain X^ m (t), taking the real part of equation 
(2.132), we obtain 
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K al K,a) I 

X im (t) = ”4~ — * cos(mwt-a T +( -l) p it/2)+ —— cos(mwt-a T ) + 


A T 


A. 


mwK, 


+ "Af 1 X m cos ("“ t + 1 + (- 1 ' 1 ~ a T > 


This gives 


X 


K ia K) J 


lm 


K l X m P 

I m sin(mwt-a T )+ — - (u> r -(-l) mu>)cos(mwt-crj.) 

(2.133) 


Similarly we can obtain mth harmonic of defined as X 2 m 

from equation (2.53) as 


X. 


K ia (-l) P 


'2m 


An 


m 


cos(raot“-a T )- 


K l 1 m 

At 




w ) sin(mwt-a T ) 


(2.134) 


It is noted from equations (2.123) and (2.134) that magnitudes 
of X^ m and X 2 m are equal and they have a phase difference of 
tc/2. 


Thus the harmonic components of the pseudo rotor currents 
X^ and X 2 are computed, Using equations (2.44) and (2.45) the 
harmonic components of the rotor current can be obtained. 

Using the above procedure, the harmonics of the pseudo- 
rotor currents have been computed for a case of particular 
motor, in Sec. 2.5. 
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2.4 VOLTAGE V/AVEFOM.i AT THE STATOR TERMINALS 

In this section, the voltage waveform, in the dq frame, 
at the stator terminals has been computed from the equation 
(1.7) s 


v dl ^ r l +L ll p ^ 1 dl + ^p 1 d2 


V ql = (r l +L ll p) + Vq2 


(2.135) 

(2.136) 


The expressions for ip and i^ 2 obtained in Sec. 2.2 have been 
used in the above equations to compute the voltage v/aveform. 

The cases of stationary and rotating rotor have been dealt vdth 
separately . 

2.4.1 For the Stationary Rotor 

For the stationary rotor, substituting for pi c10 from equa- 
tion (2,3) into equation (2.115) we have 


v 


dl " ^ r i +L 1 l p N; 1 cl 1 +M t “ ( Ip 2 pi dl + Lpp i d2^ 3 


l 2 . 


r* 


Using (a = ^ rom equation (2,18), in above, 

m2 

V cll = r l i dl“ Ma i dl + ^ L ll ~ tr ^ pi dl 
Similarly, from (2.4) and (2.136) 


^ql " r l i ql HMai q2 + ^ L ll“ I 22 ^ pi ql 


(2.137) 


(2.138) 


It can be noted, that for ideal square wave current, pi^^ and 
pi 1 are ■ j, except at the end points corresponding to 
cot = o, %/3, 2%/ 3 


y • o * 
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At t 


r osg point, and V ^ has impulse. From equation (2.117) 


ic is observed that at these points of discontinuity of i^-,,7 


1’ v dl 


S (t^o) 


"'/here 


nagni 

tude 


' 7 cll 

A 1 1 dl^ t i =0 ' > 

(2.139) 

A 1 = 

^ l h “ 4° 

(2.140) 


% 


th.ere 

is an iiioulse in i 

;he voltage waveform 


because of the leakage inductance. 

It is also observed from equation (2.137) that there a 
j jmp in V rl ^ also apart from impulse because of discontinuity of 
i dl and i^rj. If by 1 A ' we refer to the change at the point of 
consideration, from (2.137) we have 


A V dl(ti=°) “ r i^ i c n( ti=o)'-Ha/3 ida^i ~ °) (2.141) 

Substituting for^i^p (t^ = o) from (2.9) into (2.141) we have 
<*V dl ( ti=°) = A 2 ^ i dl ( ti=°) (2.142) 


where 


A 


‘-‘2 




(2.143) 


Because (aaiOp) and hence this jump is because of rotor and 
stator resistances which should bo so as there is a sudden 
change of current through these resistances. 
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Similarly , for V we have from (2.133) anc’ (2.11) 


( i . =o ) 

i ' 


V oi = A 1 Aiql (t i= o) 


\ V . ( t . =o ) = /vs A i , ( t . =o ) 

•* qi i ' — ql v i ' 


(2.1/14) 
( 2 . 145 ) 


During the time ( > o) we note that i ^ and i^ are 
constants and hence equations (2.137) and (2.133) reduce to 

V dl = r l 1 dl ” Mai cl2 for l i > 0 (2.146) 

V ql = r l i ql “ Mai q2 for > 0 (2.147) 

Thus, knowing stator anc! rotor currents expressions, the 
voltages and v ^ can be computed the impulse and jump 

magnitudes can be computed from the information of discontinui- 
ties of i ^ and i^. Table 2.3 gives the reference equation 
numbers which are are to be used to obtain voltage expression 
by use of equations (2.142), (2.144), (2.146) and (2.147). 

The table is only for interval I, II and III because expressions 
for intervals IV, V and VI can be obtained from these because of 
symmetry about (ojt = 11 ). 

2.4.2 For the Rotating Rotor' 

From equation (2.33), we have 


pi 


q2 


co r 1 c!2 


ai q2 + L ^2 ° ) r i dl " pi ql 


722 
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Substituting for 'pi 2 ' in (2.136) from this equation. 


V ql ~ r l 1 ql iA a i q 2 +to r ^ Iv5//L 22 i dl +i d 2 ^ +A l pi ql 


(2.148) 


where is given by the equation (2.140). 

Similarly, from (2.34) and (2.135) we have 


V dl ” r l 1 dl ' Mai d 2 “ Hw r^C,' 2 ‘ i ql +i q 2 ^ >+A l pi dl 


(2.149) 


From (2.148) it can be said that at the end points of the 
intervals (wt = o, %/3, 2n/3) . V ^ has an impulse and a jump 
whose expressions are given as 


6 (t i =o) 


A i A 


(2.150) 


" V ql (t i sso) = r l^ i ql“ I ' la ^ i q 2 JrL 1 W r ^ 2 * i ql + ^ i q 2 ) (2 * 151) 


Substituting for A i ^ 9 in terms of A i ^ and for Ai^ in terms 
of C\ i^j from equations (2.9) and ( 2 . 11 ) respectively, in 
(2.151) we have 


^V (t i =0) = *2 ^ 1 ql (t i = ° ) 


(2.152) 


where A 9 is ’given by equation (2.143). 


Similarly from (2.149) using (2.9) and ( 2 . 11 ) 


O (t i =o) |v dl = A l Al dl (t i =o) 

AV.,(t. = o) = Ap ^ii n (t. = o) 


(2.153) 


(2.154) 
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Thus, the behaviour of v^ and v ^ at the points of 
discontinuity of i ^ or i^ can be computed from the equations 
(2.150), (2.152), (2.153) and (2.154). 


For the remaining period of time, we note that and i ^ 

are constants, Fig. 2.2. So, for t^ > o, equation (2.148), 
(2.149) reduce to 


V 


ql = r l i ql ' Mal q2 + f '“i idl + 1 d2> 


V dl r l 1 dl " Mal d2 


ito r <L“' VV 
From equations (2.42) and (2.43) 


(2.155) 

(2.156) 


i q2 “ X 1 


Ni q 

^22 9 -*- 


(2.157) 


‘d2 ~ *2 “ T 22 *61 


(2.158) 


Substituting for and i^ from equations (2.157) and (2.158) 
in equations (2.155) and (2.156) we have 


V dl = h. bl * Ha *2 ~ “"r X 1 


(2.159) 


V ql = ^ V - Ha X 1 + M “r *2 


(2.160) 


Since X-^ and Xp expressions have been computed in Sec. 2.2.2, 
the expressions for and V ^ can be computed. 



Table 2.4 lists the reference to equations numbers of 


and X^ during intervals I, II and III which are to be 
used in equations (2.159) and (2.160) to compute and V 

The table also lists the reference to changes in currents at 
t- = o, which can be used in equations (2.150), (2.152), 
(2.153) and (2.154) to compute the magnitudes of impulses and 
jump at (t^ = o) for the waveform of and V r 


ql* 


2.5 CALCULATIONS FOR A KNOWN MACHINE PARAMETERS 

For the calculation, a slip ring induction motor 
considered in [ 12 ] has been taken. The parameters of the 
machine are given in Appendix A. 

For this machine two computer programs are developed in 
DEC System 10 computer in FORTRAN. One uses the results of 
Section 2.2 to compute the sampled rotor current waveform 
in time domain. The other program computes the various 
harmonic components of the pseudo rotor current using 
Section 2.3 approach. The listings of these programs 
are given in Appendix B. These programs have inverter 
input dc current, inverter frequency and rotor speed as 
parameters . 

Table 2.5 gives the values of currents at the start 
and end of intervals I, II and III using results of Sec. 2. 2. 2 
for rotating rotor case. It is seen from this table that the 


rotor currents at the start of eactt jinteryal i.e. at (t-=0) 
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is same as that obtained at the end of the previous interval 
at - T c )’. Thus I ( j Q and I are to be computed only 

for one interval. 

Table 2.6 gives the harmonic components of the time 
solution of pseudo rotor current obtained through first 
program. The harmonic components which are obtained by 
second program using Sec, 2,3 have also been listed. Both 
those components are seen to similar and this verifies that 
both the programs give same solution. 
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Table 2.1 

Throe phase and two phase stater currents 


Interval/ 

conducting 


o < wt < n/3 


3 phase currents 


D-Q currents 


II 


Tt, V , . , 2 % 

t wt \ 


III 

< W t < It 


'al 

'bl 

cl 


•al 


•cl 


u al 

•bl 

•cl 


“ I dc 

0 

itfc 


'dc 

"dc 


= I 
= -I 


dc 

dc 


0 


L dl 


—31/2 

i ql * “V3I/2 

Cx i^^Cwt^O) * 0 

•^ ^lCutsC) ^ 3I 
i dl ~ 0 

iql “ -V 31 _ 
AidK^/3)- 31 / 2 
A iql(ojt=rc/3)"'“' ?r3 i //2 


idl = 31/2 

iql “ ^31/2 

AidKet^/a)* 31 / 2 
& iql(a>t= 27 {/ 3 ) = V*3!/2 


IV 


% < wt < 


4it 


V 

4tc x x 5-rt 
3” i wt < T" 


VI 


bit 


< wt < 2 tc 


■al 


It 


bl 

•cl 


i.i - 


‘al 

•bl 

'cl 


"al 

'bl 

cl 


•^dc 

i dl = 31/2 

0 

i ql “ If 31/2 

^dc 

Ai dl(wt=n) “ 

A iql(wt=it) "* ^ ' 

0 

=° 

idc 

iql “ V3I 

“dc 

^^1(0)1=411/3) ** 
^ i ql(wt=4u/3) SK 

‘idc 

idl = ” 3I/2 

idc 

iql * V*3l/2 

0 

A^dl (wt=5ix/3)“‘ 
A iql(wt=5ic/3)~' 


-31/2 



Initial and final values of D-Q rotor currents for stationary motor 
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Reference for the computation of stator voltage expressions 
in the case of stationary rotor 
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Reference for the computation of voltage expressions of 
rotating rotor case 
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TABLE 2.5 
*********** 


INTERVAL 

LUMBER 


♦ ♦FOR X.1 *♦ 
INITIAL J 

VALUE 1 


FINAL 

VALUE 


tOR I0C-6.A^P, | INV 


2.96724 

* 0,13596 

*3,10320 


F -0?l2l7l HZ S . 

-km : 


FOR IDC, 6 .fHP,;|NV. FRE(!,= 20 ,HZ & 


■ 

U T 0,44359 

HI -0,92998 

FOR I»C= 8 .M? |IM. 

II I'JtfiCit, 

in • : --.3 

FOR «C= 4 .AM| ?|| »V. 

if* -?:!»!« 


-tmn : 

,B o E ?sTl§s HZ ‘ 

-0.87940 

-1.39201 

r "J!i8S» B * L 

'} .369J2 
-.1 


♦♦FOR X2** 

initial 

VALUE 

RPM-5G0 . 
-1.87004 
-3.50471 
-1.63467 

RPM=300. 

"0. 2807-8 
*1.32993 
*1.04915 

RPMssSOO. 

*0. 2147b 
•1.31291 
•1.09815 

RPM~<< 

■0.66173 

■1,71524 

1.05352 


FINAL 

VALUE 


-3,50436 
— i .64 0 6 4 
1 .87002 


*1.32837 
*1 .05037 
0.28079 


*1.31125 
*1 .09975 
0.21476 


•1.71421 

* 1.05581 

0,66172 




TABLE 2.6 
*********** 


VIA FREQ. DOMAIN 
MAGNITUDE PHASE 

INVERTER - 20 .0 HZ 

9J.|42§ 
91.190? 
90.9004 
99.5510 
90.4795 
90,3585 
90.326 ? 

’ 20,0 HZ, 

102.1206 
91.1185 
90.9465 
90,5351 
90.4923 
9 0 .3516 
90,3326 

* .20.0 HZ* 

122.7928 
91.0547 
90.3977 
90.5200 
90.5057 
90.3451 
90,338? 

~ 30,0 HE. 


3 s 6.00 FREQ OF 

It *?«2? 9.81483 

| I w 11^1 0*02650 

90.4622 O.O041O 

90,3363 0, 0</23S 

90.3017 0.OO191 

’ * 6.00 FREQ of INVERTER 

192.1182 1.33076 

lb kit | 0.02489 

0.01504 

90.5234 0.00541 

90.4746 0*00421 

90.3322 0.00230 

90.3078 0*00195 

s 6.00 FREQ OF INVERTER 

3,45320 

91,0557 0.02347 

90,9864 6.0-1586 

0:00§26 

90 . 488 1 0.00433 

1ft * ilJf '0.00226 

90.3163 0,00198 


8.00 FREQ HE INVERTER s 

•mi 1. JSIfo ! 
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CHAPTER 3 

STEADY STATE ANALYSIS OF THREE PHASE INDUCTION MOTOR 
FED BY A CURRENT SOURCE INVERTER WITH A 
NONZERO COMMUTATION TIME 


3.1 INTRODUCTION 

This chapter presents the steady state analysis of the 
three phase squirrel cage induction motor fed by a current 
source inverter with a nonzero commutation time. It has bee 
shown in [l] that the rise or fall of the line currents of 
the motor during the commutation of the inverter can be appi 
mated by a function proportional to I^ c cos w c t, where w c is 
parameter dependent on the leakage inductance of the motor a 
the value of the commutating capacitor of the Inverter. It 
also been shown in [l] that the time duration of the commuta' 
'tion, denoted by t c is given by 

t_ = Tt/2w (3.1) 

c c 

The expressions for this assumed currents waveform have boon 
given in Table 3.1. Fig. 3.1 shows the waveshapes of these 
lino currents. 

The dq model of tho induction motor , with axis attached 
the stator (Sec. l.i ), is used for the analysis. The 
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transformed dq currents, using equations (4.39), for the 
three intervals of the inverter period are also listed in 
Table 3.1 and ’I* in Table 3.1, is defined by equation (2, 
The waveforms of the dq currents are given in the figure C 

The waveform of the rotor current is obtained by the 
solution of equations (1 and the phase voltages of th 

motor arc obtained through equations (1.7). 

The analysis has been done on the similar linos as th. 
followed in Chapter 2. The difference in this case being ' 
different expressions of the line currents. Sections 3.2 * 
3.5 correspond to the Sections 2.2 to 2.5 respectively. 
Section 3.2 deals with the time domain analysis, considorir 
both the stationary and the rotating rotor cases separately 
Section 3.3 deals with the frequency domain analysis and in 
Section 3.4 the voltage expressions at the stator terminals 
have been obtained. In Section 3.5 the case of the inducti 
motor introduced in Section 2.5 has been considered. 

In this chapter, wo have an additional parameter io Q . < 
the value of w c increases, the current waveshape tends to tl- 
idenl quasi square waveshape, corresponding to Fig. 2.1. II 
has also been shown In this chapter that as w c tends towards 
infinity, the rotor currents and stator yoltage expressions 
this chapter approach those obtained in Chapter 2. 



3.2 ANALYSIS OF INDUCTION MOTOR FOR ROTOR CURRENTS THROUG1 
TIME DOMAIN 

The analysis proceeds in similar lines as done in See, 
The equation of the motor performance (1*6) is solved. The 
arc solved for three intervals (I, II, III) only as the state 
waveform is symmetric about wt = it. Each interval is consi 
red in two portions. One corresponds to the commutation po 
and the other to the remaining interval period. So the fir 
interval corresponds to the time from start of interval, t^ 
less then time of commutation, t . This analysis h as bee 
done in fl2 ]. The results are reproduced here for all the 
three intervals. It should be noted that the symbols of in. 
tial currents correspond to equation (2.2). 

3.2.1 Case of stationary rotor 

The expressions of the rotor currents, as taken from (j 
for intervals I, II and III have been given below. 

3. 7. 1.1 Interval I (o < wt < it/3) 
for o < t < t 

Hen L 

Na = No (3- 2 > 

iq 2 = (lq 0 + T 3K i) oxp(-at) -V‘3K 1 [cos(w c t)-C^ sin(w c t)] (3.3) 





3.4 

I 

qo 

~\ r 3K«K. 

= — [exp(-aTj - l/2) 

(3.5) 

K 1 

- Ma> 2 l/(a 2 +w?) L 22 

(3.6) 

C 2 

= a/w c 

(3.7) 

k 2 

® t 1+c 2 exp(at c )] exp(-aT c ) 

(3.8) 

T c 

= interval period = 71 / 30 ) 

(3.9) 


and ’a*, , C / ’ are as defined in equations (2.15) and (2.23) 
respectively. 

for t < t < T = ir/3w 

V ***'”* V 

i d2 = I do ex P(“ at ) (3.10) 

i q 2 = I q 2 exp(~a( t-t c ) ) (3.11) 

where 

,; q2 = ^qo+^l) exp(-at c ) + (3.12) 

3. 2. 1.2 Interval XI (-3 - ut - T" 1 ' e * 0 - ut 2 - 3^ 

/ 

for o < t < t 

V-> 

3K 1 3 

i d2 = (l do “ exp(~at 2 )+ 2 K JL [cos(w c t 2 )~C 2 sin(o> c t 2 )] 

(3.13) 
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V 3K X V3K, 

x q 2 ~ ( ^qo + ~ 2 ~~”) exp(-at 2 )- vjr-^ [cos(u> c t 2 ) ~C 2 sin(w c t 2 )] 

(3.14) 


wnoro 


V3Kr>K 


■qo 


1 tl _ eig>(-aT e ) ^ 


3K 2 K 1 

I do = 2 C^"~ t ex P(“ a T c ) ] 


(3.15) , 

(3.16) 


for t c < t 2 < T c> i.e.((™ + t ) < wt < |— ) 


i d 2 = I d 2 GX P(-a(t 2 - t c )) 


(3.17) 


q2 * X q2 GX P(“ a ( t 2 ”t c ) ) 


(3.18) 


where 


X d 2 = (I 


3K, 


do 


3 r 1 C 2 


2 “) oxp(~at c ) “ 


(3.19) 


f3K, f3K,C 0 

1 q 2 “ (l qo + oxp(-at c ) + — 5 — - 

3. 2. 1.3 Interval III (|& < wt < % i.e. 0 < o>t 3 < §) 


(3.20) 


for o <_ cot^ < t c i,o. < ut < (^ x - + t ) 


• 2 ti 


3 3 

i d 2 = ^ X do“ 1 K l^ exp(~at 3 ) + j Kj_[ cos(oj c t 3 ) -C 2 sin(w c t 3 )] 

(3.21) 
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i q2 = U q0 - K l> exp(-at 3 ) + T | K L [co S ( Mc t 3 )-C 2 sinC^tg) 

(3.22 

where 


I do " I %T” [oxp(-aT c )-l] (3.23 

4 


1 qo = ' f § T!" [«P(-aT c ) + 1] 

(3.24 

for t < w to < § 

C -- 3 ™ 3 


1 d2 = I d2 exp(-a(t 3 -t c )) 

(3.25! 

i q 2 = x q 2 ©xp(-a(t 3 -t c )) 

(3.26) 

where 


1 d2 = (1 do ‘ 1 K l } °*P(- a V ' 1 K 1 C 2 

(3.27) 

*q2 = ^qo - K l> “ 'I K 1 C 2 

(3.28) 


The expressions of the initial and final values of currents fo 
each intervals have been summarised in Table 3.2. It can be 
scon from the table that the expressions for the final values 
of currents in an interval are the same as the initial values 
of currents in next interval. Hence, the rotor currents are 
continuous . 
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3.2,2 Case of rotating rotor 

The expressions for quasi rotor currents because of th< 
stator waveform of Fig. 3.1, have reproduced in this section 
from [l2]. Actual rotor currents can be computed from these 
quasi rotor currents by using equations (2.44) and (2.45). TJ 
symbols used below correspond to equation (2.58). 

3. 2. 2.1 Interval I (o < ut < ^0 

for o < t < t 
~ c 

Xj(t) = exp(-at) [A q1 sin(m r t)-fB ol cos(u> r t)] - K 4 + 

4- K^a[m c sin(w c t)+K 2 cos(y c t) ]+K^w [co cosa> c t-K 0 sinw 

(3.29) 


whore 


B ol 

— V * /-* 

~ yv lC ' r U 1 

(3.30) 

C 1 

= K 4 - iqK^a - KjA 

( 3 .31) 

A 1 

01 

= (X 10 +aB ol -K 5 a W ®+K 5 K 2 A)/ Ur 

(3.32) 

2 

= (w 2 -i-a 2 ~w 2 )/2a 

(3.33) 

K 3 

= l/(w 2 -ei<|) 

| 

(3.34) 

K 4 

= ][ § K X I (r3w r a)/(a 2 +m 2 ) 

(3.35a; 
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K, 


M a /L 22 


(3.35b) 


K 5 = f 3K 1 IK 3 /2a 

X 10 = w r X 20~ aX 10 + K ll 

and 


(3.36) 

(3.37) 


X 3 (t) sxp(-at) [A q 2 sin(w r t)+B Q 2 cos ( w r t )3 + K ( 

K 5 w r (a) c sin ‘V + K 2 cos t )3 


whers 


(3.38) 


3 o2 = X 20 + C 3 

C 3 = KgKgWr-Kb 

A o 2 = (X 2O +aB o2 +K 5 w r a) c )/t0 r 


(3.39) 

(3.40) 

(3.41) 


K 6 = fSI^lCw -f 3a )/2( a 2 +o) 2 ) 

**• X* 


X 20 ““r X 10 ~ aX 20 “ 2 lx l 


K, I 


For X 20 and x io expression is given by 


! X 20 j ! l/2-(Y 4 Y 2 - Y 1 Y 3 ) - f| +(y 


'42 I 1 I 3 / ” 2 


! x io | j V l -(YaWd 


5’ -C Y 4 Y 2- Y 1 Y 3>! 


(3.42) 


(3.43) 


3 Y 2 +Y lVi ! Z. 


j 2. 


(3.44) 
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where 


Z 1 ~ C 2 Y l +C l Y 2" K 4 +K 5 aco c ~ K 5 a) c K 2 

(3.45) 

Y i = oxp(-at ) sin w t 
x c re 

(3.46) 

Y ? = exp( -at) cos w t 
c re 

(3.47) 

Y 3 “ °xp("a(T c -t c )) sin[ Wr (T c -t c )] 

(3.48) 

Y 4 = ex -°[ _3 ( T c “'t c )] cos [s> ( T -t )] 

x c o 

(3.49) 

C 2 = K 1 I+aC i"~K 1 a'o2 + K 5 K 2 w2] /w r 

(3.50) 

Z 2 = C 4 Y l +C 3 Y 2 +K 6“ K 5 w r C 2 

(3.51) 

C 4 = [ c 3 a “ 1 K l I+K 5 0J r w c 1 / co r 

(3.52) 

Z 3 = C 5 Y 3 +(Z 1 +K 4 )Y 4~ K 4 

(3.53) 

C 5 = [~ K 1 I+K 4 a+w r Z 2 ]/co r 

(3.54) 

Z 4 = 

(3.55) 

C 6 = ^“ Z i w r“ 1 K l I “ K 6 a ^/° J r 

(3.56) 
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The order of calculations of the constants for calculation of 
^ 10 ono X 2 q Is to to Y^j to '^“2.* 

Z„ and Z A . 

for t c < t < | 

:< x (t) = cxp[-a(t-t c )] [C Ql sin(w r (t-t c ) + D q1 cos(o> r ( t-t c ) ]-K 

(3.57 

X 2 (t) = exp[-a(t-t c )] [0 o 2 sin(ui r ( t-t c ) )+D q2 cos(w r ( t-t c ) )+K^ 

(3.58’, 


where 



D ol 

ii 

X 

t— 2 ** 

+ 

(3.59) 

C ol 

= [Xi + aD ol ]/u r 

(3.60) 

x i 

" X 20 Y 1 +X X0 Y 2 +Z 1 

(3.61) 

a 

X i 

= u r Xt, - aXi - V| K X I 

(3.62) 

X 2 

= _X 10 Y 1 +X 20 Y 2 +Z 2 

(3.63) 

D o2 

h 

! 

X 

(3.64) 

C o2 

= [X^ + aD o 2 "j/ w r 

(3.65) 

X 2 

“ -“r X i - oX 2 - 1 V 

(3.66) 
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3. 2. 2. 2 Interval II (i.e. (i.e. | < t < ^ or o <_ wt 2 < 
for o < t 0 < t 

— 2 C 

X JL (t) = exp(-at 2 ) [A q3 sin(u> r t 2 )+B o3 cos(o3 r t 2 )] + 

~K 7 +Kq[w c sin(w c t 2 )+K2 cos(w c t 2 )] + 

+ K 9 [o) c cos(w c t 2 )-K 2 sin(to t 2 )] (3.67) 


X 2 (t) = exp(~at 2 ) [A o4 sin(o) r t 2 )+B o4 cos(a> r t 2 )] + 

K^o""Xn[a>c sin(0) c t ? )-rk 2 cos(co^t 2 ^) + 

- k 19 [w c cos(u) c t 2 ) - K 2 sin(co c t 2 )] (3.63) 


where 





B o2 

= X 10 + K 7 - *8*2 " K 9 W c 



(3.69) 

CO 

o 

< 

= [X 10 +aB o3 »K 8 w^rK 9 K 2 u) c ]/o.> r 



(3.70) 

X 10 

= w r X 20 - aX iO - K 1 I 



(3.71) 

K? 

= (f3IK 1 a)/(a 2 +u>^) 



(3.72) 

K 8 

= [(V3I K ia /2)-(3I K x w r /2)] 

(K 3 /2a) 


(3.73) 

k 9 

= f3I K 1 oj c K 3 /4 


\ 

(3.74) 
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B o4 

^O^lO^l 1 K 2 " K 12 u c 

(3.75) 

A o4 

[^20 +aB o4 +K ll u) c“ K 12 K 2 w c^ W r 

(3.76) 

*20 

3 

** a ^20 ** J 

(3.77) 

O 

1 — i 

= V3I K,w^/(a 2 +ca^) 

1 1 r 

(3.78) 

K n 

= C ( 3K x Ia/2)+(f3K 1 w r l/2)](K 3 /2a) 

(3.79) 

K 12 

= 31 K x w„ K 3 /4 a 

(3.80) 


The values of X^q and X 20 is determined from matrix equation 


■20 

! 

1 

1 r 1 

i ! 2 - 

[ 

(Y 4 Y 2 -Y 1 Y 3 ) 

- f i + ( Y 3 Y 2 +Y 4 Y d | ’ 

00 

N 

1 

J 

I 

, f3 
j 2 “ 

C Y 3 Y 2 +Y 4 Y 1> 

5 - (Y^-YjYg) i 

( 

i 

j Z 7 


(3.81) 


The constants are to be solved in following order to Kg, 
K.y to s to C 7 to Cj^q s Z^ f Zg y ^ 11 * C<i 2 * Z 7 and Zg 
where 

C ? = K ? ~ KqIC, - K g w c (3.82) 

C 8 = [ aC 7 -(Y'3K 1 l/2)-K 8 u^+K 9 K 2 u) c 3/w r 


(3.83) 
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C 9 = ~ K 10 +K 11 K 2 +K 12 0) c (3.84) 

C 10 = f'(” 3K 1 I / 2 ) + aC 9 +K 11 a)^-K l2 K 2 w c ]/oo r (3.85) 

Z 5 = C 8 Y 1 +C 7 Y 2 “ K 7 +K 8 W c ~ K 9 K 12 (3.86) 

Z 6 = C 10 Y 1 +C 9 Y 2 +K 10“ K 11 w c +K 12 K 2 (3.87) 

C ll = ' LZ 6 w r ~ V"3 K 1 I-hK 7 a]aJ r (3.88) 

C 12 = ^ Z 5 w r“ K 10 a ^ w r (3.89) 

Z 7 = C 11 Y 3 + ^ Z d +K 7^ Y 4“” K 7 (3.90) 

z 8 = C 12 Y 3 +(Z 6“ K 10 )Y 4 +K 10 (3.91) 

For t c £ wt 2 £ § 

/l(t 2 ) = exp(-a( t 2 -t c ) ) [C o3 sin(w r (w r ( t 2 -t c ) ) + 

+ D o3 cos(u r (t 2 -t c ))]-K 7 (3.92) 

~ ex PC-a(t 2 - 1 ; c ) ) [C o4 sin “ r (t 2 -t c )) + 

+ D q4 cos(ts r (t 2 -t c ))] + K 1q (3.93) 

where 

D 03 = 4 + k 7 


(3.94) 
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n 

°o3 

= [ x i + aD o 3 l 

(3.95; 

D o4 

= 4 - K 10 

(3.96) 

C o4 

= [4 + aD o4 ]/w r 

(3.97) 


'A\ and can be obtained from equations (3.67) and (3.68) by 

• « 

substituting t^ = t c . and XT, can be obtained from equatio 

(2.50) and (2.51). 

3.2 .2 .3 Interval XII (i.e. u>t b te or o js u>t 3 S, 

for o < t < t 

X 1 ( t) = exp(-at c ) [A o5 sin(o> r t 3 )+B o5 cos(u> r t 3 )] + 

+ %3'" K i4[ w c sin(to c t 3 )+K2 cos(co c t 3 )] + 

“ K 15 Ew c cos(o) c t 3 ) - I <2 sin(t«J c t 3 )] (3.98) 

~ exp(-at c ) [A o6 sin(o) r t 3 )+B Q6 cos(a> r t 3 )] + 

+ K 16 +K 17^ u c sin ( (U c t 3)' ! '% cos(w c t 3 )] + 

- V 3 K 15 [w c cos(w c t 3 ) + sin(oj c t 3 )] (3.99) 

where 

B o5 = X 10 ™ K 13 + K 14 I<C 2 + K 15 


(3.100) 



3.15 

A o5 = ^10 +aB o5 + K 14 a) c“ K 15 K 2 u) c^ w r (3.101) 

^10 = w r"^20 ~ a ^io “ K^f31 (3.102) 

iX j_3 = [V*3K^l/2j (V*3to r "a) /( (3.103) 

K 14 = ECV 3 ^ 1 a/2) +(3K.,Iu> r /2)] (K 3 /2a) (3.104) 

K 15 = f (3.105) 

B o6 “ X 20“ K 16“ K 17 I< 2 + '^ 3K 15 w c (3.106) 

A oo = ^0 +a3 o6~ K 17 W c~^ 3K 15 K 2 Cl> c ;i/w r (3.107) 

* 

X 20 = -“r X 10 " aX 20 ( 3.108) 

K 16 “ (V’3K 2 I/2)(iii r +ir3a)/(a 2 'i-u 2 ) (3.109) 

K 17 = [(f3K lMr I/2)-(3K 1 aI/2)](K 3 /2a) (3.111) 

For the values of and X-jq, we have 


X 20 


f 

2 “ ( Y 4 Y 2“ Y 1 Y 3^ ^2 + ( Y 3 Y 2 +Y 4 Y 1^ 


j 

z 12 

X 10 


T | ~ (Y 3 Y 2 +Y 4 Y 1> k “ ‘ Y 4 Y 2- Y 1 Y 3’ 

1 

1 

I 

Ni 

l 


(3.111) 
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j. he constants are to be evaluated in the following order 


K 1 t0 K 3> K 13 t0 K 17> Y 1 t0 Y 4- C 13 t0 C 16- 
^11 and Z-,,-,, where 


Z 9Sf Z 10 ,C 17> C 18 , 


n 

^7 o 

JL 

= -*K 13 +K 14 K 2 + K 15 a) c 


(3.112) 

C 14 

= [ -f3K 1 I+aC l3 +K 14 w^-.K 15 K 2 a) c ] /o> r 


(3.113) 

C 15 

= ~ K 16“ K 17 IC 2 + > r3K 15 co c 


(3.114) 

C 16 

= CaC 15 -K 17 u|-V3K 15 K 2 u c ]/a) r 

- 

(3.115) 

Z 9 

C 14 Y 1 +C 1 3 Y 2 +K 1 3 " K 1 4 w c +K 15 K 2 


(3.116) 

Z 10 

= C 16 Y 1' < ' C 15 Y 2 +K l6 +K 17 a> c + ^ 3K 15 K 2 


(3.117) 

C 17 

~ f ^*10^r” (Y"3K. x l/2 ) “K^ga 3 /to.^ 


(3.118) 

C 13 

= [“Z< 5 w r +(3K 1 l/2)-K 16 a]/w r 


(3.119) 

Z 11 

= C 17 Y 3 +Y 4 (Z 9 -K 13 ) + K 1 3 


(3.120) 

Z 12 

= c 18 y 3 +(z 10 -k 16 )y 4 +k 16 


(3.121) 

For 

t c i < | 




X x (t) = exp(-a(t 3 -t c ))[C 05 sin(o> r (t 3 -t c )+D o5 « r (t 3 -t c ) )]^ 13 


(3.122) 
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X 2 (t) = exp(-a(t 3 ~t c )[C o6 sin(w r (t 3 -t c ))+D o6 cos w r ( t 3 “ t c ) )-3 +K 

(3.123 


where 

°05 ’ X i - K 13 

i -'o5 ~ ^1 + a ^o5^ w r 

D o6 = *2 “ K 16 

C o6 = [ *2 + aD o6 ]/t0 r 


(3.124 

(3.125 

(3.126 

(3.127 


The expressions for and X£ can be obtained from equations 

* 

(3.98) and (3.99) respectively by substituting t 3 = t c . and 
X4, can be obtained from equations (2.50) and (2.51) by sub- 
stituting for X^ and X^. A computer program has been develope- 
to calculate initial and final values of currents for these 
intervals in the rotating rotor case and the boundary condi- 
tions are found to be matching (Sec. 3.5), 

3.2.3 Case of u> — v 00 

V 

As a) increases, the stator current waveform of this 
chapter, (Fig. 3.1), tends towards the waveform considered 
in Chapter 2, (Fig. 2.1). In this section it has been shown 
that as w c tends towards infinity, the expressions for rotor 
currents obtained in Sec. 3.2.1 and 3.2.2 reduce to that 
obtained in Secs. 2.2.1 and 2.2.2. This reduction is shown 
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for interval I for both tho sections. The results for other 
intervals can bo shown by proceeding on similar lines. 

3. 2. 3.1 Stationary rotor case 

As u) <» wg note from equation (3.1) 

V 

t = 0 (3.128) 

So tho interval o < t < t of Sec. 3. 2, 1.1 reduce to a point 
t = o and the period t < t < T to tho remaining period of 
interval I except t = o s as is the case in Sec. 2. 2. 1.1. 

Also, as m co 
c 


K = ML- = K* [From (3.6), (2.13)] (3.129) 

x l 2 2 x 

C 2 = 0 [From (3.7)] (3.130) 


K 2 = exp(-aT c ) = KJ> [From (3.8) ,(2.23) ,(3.130)] 

(3.131) 


Thus from (3.129) and (3.131) I^ 0 in equation (3.4) reduce to 

that in equation (2.21) and I of equation (3.5) to that in 

qo 

equation (2.22) . Also I^ 2 of equation (3.12) reduce to 
Iq2 “ (lq 0 + V^K^) and this equation is similar to (2.12). 

The expressions of currents in Sec. 3. 2. 1.1 become 


i d2 = tio ox P<-at) 
V “ r q2 0 X P(*at) 


[From eqn. (3.10)] 


(3.132) 


[From eqn. (3.11)] 


(3.133) 
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Those expressions are equivalent to equations (2.16) and 

(2.17) os P 2 of (3.133) is already shown to reduce to that in 

(2.12) and I^ 2 1^ 0 f r om oqn. (2.10) in equation (2.16). Thus 

the expressions obtained in 3. 2. 1.1 reduce to that obtained in 

* 

2. 2. 1.1 as to r approaches infinity. 

3. 2. 3. 2 Rotating rotor case 

Here it will be shorn that for co approaching infinity, 

V 

the rotor currents expressions of section 3.2,2, 1 reduce to 
that obtained in Section 2. 2. 2,1. It has already boon shown 
in 3. 2. 3.1 that sub regions of the interval taken in 3.2.2. 1 
reduce to that taken in 2. 2. 2.1 as w c --- °°» 

For t = o 

So, eqn. (3.29) is for a point t = o. Thus from (3.29) 
X x (t=o) - B q1 - K 4 + K^aK 2 + K g uF; (3.134) 

Substituting for B oX from (3.30) in (3.134) 

X x (t=o) = x i 0 +C i” K x' l_K 5 K 2 a + (3.135) 

Substituting for C x from (3.31) in (3.135) 

X 1 (t=o)* = X 10 (3,136) 

as should bo 

Similarly eqn. (3.38) reduce to 
X 2 (t=o) = * 
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For X 1Q and X 2Q 


As m --*■ «>, since t = 0 , we have 
Y^ = O [from (3.46)] 

Y 2 = 1 [from (3.47)] 

Y ^ = exp(-al c ) sin u r T c [from (3.48)] 

Y ^ = exp(-aT c ) cos to r T c [from (3.49) j 

C x = K 4 [from (3.31) , (3.33), (3.34), (3.36)] 

Z 1 = 0 [from (3. 45), (3. 36) ,(3.34)] 

c 2 = [- ^I+al^] /u r [from (3.50)] 

C 3 = [from (3.40)] 

C 4 » [**K 6 a - K 1 l]/w r [from (3,52)] 

Z 2 = 0 [from (3.51) , (3,136) to (3.144)] 

= [- K 1 I+K 4 a]/m r [from (3.55)] 

C 6 “ ^ I ^i 1 " K^a]/^ [from (3.56)] 

Z 3 = C 5 Y 3 + K 4 Y 4 - K 4 [from ( 3. 53) , ( 3 .141) ] 

Z 4 “ Y 3 C 6“ Y 4 K 6 +K 6 [from (3.56) , (3.145) ] 


(3.136) 

(3.137) 

(3.138) 

(3.139) 

(3.140) 

(3.141) 

(3.142) 

(3.143) 

(3.144) 

(3.145) 

(3.146) 

(3.147) 

(3.148) 

(3.149) 
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and from (3.44) using (3.136) and (3.137) 

(3.15C) 

Thus it is observed that and of eqn . (3.146) and eqn. 
( 3 . 147 ) correspond to and obtained in equations (2.72) 
and (2.75). Also \' 3 and Y 4 of Sec. 3. 2. 2.1 in reduced form 
for (*) -■» °° in equations (3.138) and (3.139) correspond to Y^ 
and Y 4 of Section 2. 2. 2.1, i.e., equations (2.68) and (2.69). 
Thus Z 3 and Z 4 from (3.148) and (3.149) are seen equivalent to 
Z 3 and Z 4 respectively, as in equations (2.71) and (2.74). 

Thus equation (3.150) give the same value of X^q and 
give the same value of "^20 an< ^ ^10 as obtained from equation 
(2.76), implying the similar initial conditions. 

For t > o 

The equations (3.57) and (3.58), using (3.123) can be 
seen to reduce as 

X-,(t) exp(-at) [C pi sin V + D ol COS (i) ^ b 3 ' ( o * 13 1 ) 

X 2 (t) - exp(-at) [C o2 sin a> r t + D o? cos a> r t]+K 6 (3.152) 

It can also be shown from equations (3.59) to (3.66) that as 

0} " * co 

c 
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X i - X 10 

(3.153) 

“2 a 20 

(3.154) 

X i = “r X 20‘' a/v 10 " ’I K 1 I 

(3.155) 

X 2 = -tlJ r X 10 -aX 20 - 1 V 

(3.156) 

U ol = X 10 + K 4 

(3.157) 

C ol= + aD od/“r 

(3.158) 

D o2~ ^0 " K 6 

(3.159) 

C o2 = [% + aD o2 ]/a, r 

(3.160) 


It is clear from equations (2,59), (2,60), (2.61), (2.66) and 
the above equations that 

t,2= D o2 < 3 .161) 


and 

C o2 = C o2 


(3.162) 


Thus equations (3.151) and (3.152) correspond to equations 
(2.56) and (2.62). Thus, as co ~r 08 , the results obtained in 

C 

Sac. 3, 2. 2.1 reduce to that obtained in Sec. 2. 2. 2.1, as it 
should be. 
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3.3 ANALYSIS OF INDUCTION MOTOR FOR ROTOR CURRENTS THROUGH 
FREQUENCY DOMAIN 


The analysis to obtain the rotor currents through frequenc' 
domains has two steps, as is done in Soc. 2.3. First is the 
calculation of the harmonic components of tho stator current. 
Thon the harmonics of quasi rotor current, and ^ are com- 
puted by tho use of equations (2.52) and (2.53). The harmonics 
of actual rotor currents can be obtained from equations (2„42) 
and (2.43) after wo have computed the harmonics of stator and 
quasi rotor currants. 



For the harmonic components of the stator current, 
and i ^ can be written in Fourier series as 

CO oo 


Fig .3.1, 




cos nwt + 


E 

n=l 


d„ sin ntot 
n d 


(3.163) 


i 1 = £ C cos nwt + 

q n-1 q 


OO 


L 

n=l 



sin nwt 


(3.164) 


From Fig. 3.1, it can be said that for even values of n, the 

constants C , C , d„ and d are zero because of mirror image 
n d n q n d n q 

symmetry. For odd, values of n, 



[-sin 



2_ 

nw 


1 

(nw+m c I 


1 

(nw-w c ) 


+ 


, 1 ' 1 
+ 2 Tnw^wTJ 

V 1 





)+cos ( 


rotw 



+ § 5 ); ] 

(3.165) 
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°n d - T 1 TS5^T C T " XreJHiJX x Sj - n( 23; 


T~ } T 


, ^ „ / mtu> , 2nn % ’ n 
h Sml!T“ -r “ ?s™j J 

*v^ 


~ - 2 nr: . 

C n ~ yi s 3 d n, 
q s c 


2 


nn 


*n = T3 sin ‘ 3 V 


(3.166) 

(3.167) 

(3.168) 


It can be seen from equations (3.166) to (3.168) that for n to 
be multiple of 3, these constants are zero. Thus from now* 
when we refer to harmonics we mean onl y odd and non tr ipi ents. 

Thus* writing the mth harmonic component of i^ and i . 
as 


i ila = I ra, costrocot + p ) 

Q G 

d qlm ~ cos ( rnw 't + P m q) 


We note from equations (3.165) to (3.168) that 


jl« j 

= jl* 

[ = 

1 1 » j 




n d 

m q 






p m = 
m q 

(P_ - 

m d . 

TC\ 

? ) 

for 

m 

= 1,7,13, 

19 ... 


+ 

1t\ 

2' 

for 

ra 

= 5,11,17 

, 23 , ... 

q 

d 





If we 

define 

•p* 

as 

in 

equation 

(2.123) then 


(3. 169) 

(3.170) 


(3.171) 

(3.172) 


(3.173) 
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“dim m 


= 1^ cos(mmt + p m ) 


(3.174) 


qlff. = l m c ° s < nMt + ?m d +( - )P 5> 


(3.175) 


where 


^ +<£) 


2 xl/2 


(3.176) 


= tan^ (-d m /C ) 


(3.177) 


It should be remembered here, that m refers only to odd and non 


riplent value. Solving (2.52) and (2.53) for mth harmonic of 


X-, and X 2 respectively, as done in Sec. 2.3 we can again visua- 


lise the picture as in Fig. 2.3. For this case we have 


K^-DP 


II sin(mwt + - a ) + 

Ji Ai.1^ X 


+ —i— ® ( 0 ) —(— 1)P raw) cos(raot + p 
a T •*- m d 


(3.178) 


K,al 


1 m 


sin(mwt + - ctj) + 


K.I» ‘ 

+ — ™J3 (co -(-1)P ma>) sin(mwt4-6 -a-,) 

x m I 


(3.179) 


where and acj are given by equations (2.129) and (2.130). It 
can bo noted from equations (3.178) and (3.179) that the 


magnitude of harmonics X^ m and X 2 m are equal and they have a 
phase difference of tc/2. In Sec. 3.5, harmonics of X^ and 


have been computed for various values of w and w T for a case ol 

V * 


induction motor. 
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Case of lo — *■ <*> 

As ti> c increases, the stator current waveform of this 
chap car (Fig. 3.1) tends towards the waveform considered in 
Chapter 2, (Fig. 2.1). So, in the limit of w tending towards 
infinity, the result of this section tond to that obtained 
in Soc. 2.3, as is shown bo low 



(3.1SO) 

(3.181) 


1 * \ J * 

Lt C = a (3.132) 

O - co n d n d 

c 

Thus, from equations (2.126), (3.176), (3.180) and (3.182) 


Lt 

iJ c 





from equations (3.177) and (3.180) 


(3.183) 


,3 r , = 0 (3.184) 

u *d 

Substituting for IJ and p from equations (3.183) and (3.184) 

rn lfl d 

in (3.178) and (3.179) we see that those expressions of har- 
monics of quasi rotor currents reduce to that obtained in 
Soc. 2.3, equations (2.127) and (2.131). 
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3.4 VOLTAGE WAVEFORM AT THE STATOR TERMINALS 


The expressions for the voltage waveform at the stator 
terminals for the stator dq currents as in (Fig. 3.2), can be 
obtained by proceeding similar to Sec. 2.4. For this analysis 
each interval is divided into three regions. They correspond 
to time (t- < t_), (t. = t_) and the remaining period of the 
interval. In this section the cases of the stationary and the 
rotating rotor have been dealt with separately. 

3.4.1 For the stationary rotor 


The equations for stators voltages for this case of 
stationary rotor have been computed in Sec. 2.4.1. They are 
equations (2.137) and (2.138), i.e.. 


V dl = r l 1 dl ~ Ma i d2 + ( L ll " Lj 2 ^ pi dl 

M 2 

V ql = r l i ql - Ma i q2 + (L 11 “ E^ pi ql 


(3.185) 

(3.186) 


Since i^ and i ^ are continuous, (Fig. 3.2), there is no 
impulse in this case. 


At t. = 0 



The jump in V^ and V at t^ =0, given by 
can be computed from (3.185) and (3.186) 


A'V^i and 
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'dl 


= (L 


11 


) 

h2 


U( 




‘dl y 


t.=o 



(3.187) 


(3.188) 


Table 3.1 can be tjsed to compute the changes in slopes at 
= o. From Table 3.1 it can be shown 

^eWk=o * 0 < 3 ' 1S9 > 

i i 


' (pi 



V 0 


0 


(3.190) 


where 1 = 1,2,3, ... 


and so V dl and V ' ^ continuous at t^ =0. 

For o < t. < t 
1 _c 

The equations (3*185) and (3,186) can be used to compute 
the expressions of the voltages. Table 3.4 summarises the 
references to the expressions of -dl’iql’-Ms ^q2 are 

required for the computation of voltages. 

At t. = t 
2 . c 

» -cuscrnr r +A*?ar - t 

It is seen from (Fig. 3.2) that there is sudden chance of 
slopes of i^ and i •, at this point. Hence, there is a jump 
in the voltage expressions at t - = t , as can bo seen from 
equations (3.137) and (3. 188). 
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From Table 3.1, using t_ = from equation (3.1), wo have 

AJJj _ 

C 


dl^ j 

tt=t c ) “ 0 

(3.191) 

ol> 

<t«t e ) =^ 3I “c 

(3.192) 

dl>j 

( t 2 =t c )“ " 1“ “c 

(3.193) 

) i 
ql'j 

!(t,=t )“ ,J, c 

(3.194) 

dl^ 

_ 31 

( t 3 =t c ) 2 c 

(3.195) 

ql^ 

i 

i f31 /% 

|(t 3 *t c } 8 “ ‘sr w c 

(3.196) 


Equations (3.19l)to (3.196) can bo used in equations (3.137) 
and (3.133) to compute the jumps in and V ^ at = t c » 


for t < t. < T„ 
c i c 


During this interval i H ^ and i r( ^ are constant. Hence, 
equations (3.185) and (3. 136) reduce to 


v ai = - Ma ta 


V ql " r l i ql " Ma 


(3.197) 
(3. 193) 


Tabic 3.3 gives the refer c-nces for the computation of the 
voltages . 
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3.4.2 For the rotating rotor 


From equations (2.148) and (2.149), using equations 
(2.42) p (2.43) to elliminate for i^ and we have 


V di — ” M 3.Aj^ — Mu> r X^ + A^pi^ (3.199) 

V ql = ‘Vql " + '*>r X 2 + V 1 ,! <3.200) 


where A^ and A 0 are defined in equations (2.140) and (2.143) 
respectively. 


Since i^ and i , are continuous. Fig, 3.2, 
are also continuous, equation (2.45) and (2,46), 
impulse in the expressions of and Vg^. 


and X 1 ,X 2 
there are no 


At t ± - 0 


It has been seen in Section 3.4.1 that at t^ — 0, the 


changes in and ’pi^’ are zero and hence and V 


qi 


is continuous at t. = o. 


For oft. < t 
__ i c 


The expressions for X^ and obtained in Sec. 3.2.2 to- 
gether with Table 3.1 is used in equations (3.199) and (3.200) 
to obtain the expressions for voltages. Table 3.5 lists the 
references to equation numbers for and 
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Fo- r t ± = t c 

From equations (2.199) and (2.200) because an 

X g are continuous 


'^dl^i~^c^ “ A^tpi^i) 


l u dl' 


'■ V ql (t i =t c) = 


t.=t 
1 c 




(3.201) 

(3.202) 


The e quations (3.191) to (3.196) can be used in equations 
(3.201) and 9(3.202) to obtain these discontinuity of voltage 


For t < t, < T 
c 1 c 


In this region i^ and i ^ are constant and hence equa* 
tions (2.199) and (2.200) reduce to 


V dl = *2 *dl ~ MaX 2 - Mw r X l (3.203) 

V ql = A 2 ^1 “ MaX l + Mw r X 2 (3.204) 

Table 3.4 lists the references to the equations which give th< 
expressions for X-, and iq]. during their region, and 

can be used to compute v^ and v . . 

3.5 CALCULATIONS FOR A KNOWN MACHINE PARAMETERS 

Using the results of Sec. 3.2 and Sec. 3.3 two computer 
programs are developed to compute rotor currents. The listing 
are given in Appendix B. 
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Using those programs the data for Tables 3.5 and 3.6 arte 

obtained. It is seen from Table 3.5 that the I , and I „ 

• do qo 

values need to computed only for one interval. Table 3*6 
shows that the solutions obtained through both the programs 
using Sec. 3.3 and Sec. 3.4 results are similar. It is also 
seen from Table 3.6 that as w increases, the results approach 
to that obtained in Table 2.6 for ideal inverter current. 




Three- phase and two phase stator currents 
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Initial and final values of D-Q rotor currents f 


I „C\! 

= tr 


u 

H 

4“ 

r | 

i 

O 

? 

u 

c 


H 

H 

Cl 

CJ 

* 

n: 

X 

if 

s 

O 

Cl 

CL 

rH 

X 

X 

i 1 

r j 

o 


I *i . 

I 

CO 


4V 


k fMi *3 

fx. p 

CO | 


I— if 

^ ! 

i , p i ^ 

™ p 


:iv 

CO 1cm 


r-T 

3 

Osjf 'sj- 
U£ JO 
CO JCN 

I 


! 

! £!° 


04 i 'T X 

^ Q O 

a • 


£ 

O k| o 


£* & 

| CO 5 cm 


M? 


>HS 

& 

<*> !cm 


CO 

\ 

^ M 

V / S *-4 

T* 

■H > 

3 cc 

Vf f-t 

~ ‘Z. 

° t— 1 


ro 

M 


h~i 


t — f 


M 

CM 

vj 

*4 

< 

Mt 


*P 

VI 

&' 

W 

3 

\ 

•> 

cC 

JTJ 

CO 

\ 

H 

fc: 

£ 

*—* 

W 

£ 

CM 

M 



C) 

C 

o 

P 

CO 

<0 

o 

Q, 

X 

o 

O 

£r> 

fC 

-p 

P 

O 

> 

M 

O 

P 

03 

p 

CO 

p 

o 

CO 

* C 

CO o 

P 

P O 
0 CO 07 

ip P f0 

P 3 u 

ns o. 

^ i g ^ 

o o 
o P 
o 
O M 
p 

p > 

p 

U (0 

o c 

P O 
P 
O P 
u ro 
C P 
Q 07 
P 

O CO 
P 

OP 

cC o 


5 r~i 


p o? 

v ‘ — 3 * p 

Cl * CU * 

tro ctco 

O ' •' (*J V „ 


NO 

• CM 

c * 

CTCO 

O 


14 ^ 1 ? 


; <t> . p 


o 

« ^_j 

a • 

G'oo 


o 

• I — i 

c - 

troo 
o — 


lO 

• CO 

c * 

vTCO 

O — 


o 

p 

p p 

a • 

P CO 


X> rH 
05 * 

P CO 


o 

p 

_c p 

o • 

P CO 


Oi 

O'* 

« p 

C * 

croo 

O — 


ON 

* P 

c - 

trco 

O — 


\D 

O'* 

• p 

c • 

creo 


CO P 

H a 


P 

o 

o p 

c * 

CTCO 

O — 


CO 

o 

o 

\C * 

creo 

o— 


lO 

o 

• p 

c * 

cr co 
o — 


p 


LJ cm 

* ni J-' 


O 


5 °! 


J4 


Kfcj ° 


i jiOj <h 

r r 

| *H| •! *H 

+> !C0 -o tr 

I - t — 'f >rJ ‘ 

V] . 

L <C] u x; 

1° jew O C 
lajf tu ra 

o 

j ■ J & i w 
Sjc U O C 1 O 

|> C -P ^ o O o 

j£j 9 'tj 'p c 
j- -H O X W O O 

f£ g 1 > O <o n H 

(H 


1 


i • w 

j C • 

{ treo 
o — 

1* 


• (N 

c * 

creo 

o— 


N 


o f 

-P-4 

C cs 

P Ji 


o 

p 

p p 

05 • 

P 00 


M 


» p 

C * 

creo 
o — 


a p 

C • 

creo 

o— 


o 

p 

p p 

CO * 
P 00 


OJ 

•Csj 

a • 

O'* 00 

o — 


• C"4 

jTQ « 

treo 

o— 


o 

p 

p p 

c • 

P 00 


ip 


M 

Hi 


Table 3.4 

RofereiiGo for the computation of stator voltage oxorossions 
of rotating rotor case 
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TABLE 3,5 


INTERVAL 

HUMBER 


I 

II 


** FUR XI ** 
INITIAL FINAL 

value value 


** FOR 
INITIAL 
VALUE 


X2** 

FINAL 

value 


FOR lDC«t»,AMPjINV. FREQ ,=20 . HZ; RPM*500. ? WC=500 , 

I 3,26250 0,69f§4 -1 ,09295 -3.37000 

II 0.66468 -2.57277 -3,37165 -2.28387 

HI -2.57730 -3,26243 -2.27890 1.09291 


II 0.92027 -0.65605 

III >V> -1,77973 

FOR I DC *4 .AMP I IN V. FREQ. =15 


1. 07192 0.41953 

0.41779 -0,65515 

III -0,65713 -1,07490 

FOR IDC*6,AMPj INV. FREQ. =20 


,HZ;RPM=600,?WC=860, 
0.935 $.2 -1.52152 

-1.52376 -1,56027 

-1.55888 -0.03512 

,HZ;FPM=200.?<*C=/50, 

-0.13816 -0.99872 

-0.99998 -0.86266 

-0.86181 0,13816 

.HZ?FPM=4OO.?4C=70O. 


= I 
III 


1,94711 
6.74558 
-1 .20152 


0.74907 

■1,19805 

-1.94707 


-0,26321 

-1.81783 

-1.55462 


-1.81587 
-1,55661 
' 0.26320 



TABLE 3.6 

* $ t * '(t { £ $ $ 4: * 



FREQ. 330 . HZ & WC= 800 
1.34655 85.647^4 

f. 01983 22.87026 

0.01 t 07 -4.92045 

0.00289 -o2.9u4b6 

0.00182 266*05194 

0.00066 193.78426 

0,00035 148.23648 

FR?Q, =f5.MZ i W* 1O00 
1.00010 99.05103 

0.02200 64.48393 

0*01285 53.38765 

0.00448 31 ,0y987 


w * o.«wri a 1. » 

85.64450 
22.83690 
-4.90470 
-62.91205 
26o. 12643 
194.26089 
148.15403 

dc s 4. amp; in v 

99.04695 

£ it " A £■ m ti 


FREQ, 


—————FOP XI— — — — — 

VIA TIME DOMAIN VIA FREQ. DOMAIN 

MAGNITUDE PHASE MAGNITUDE PHASE 


FOR RPM = 100. ? IOC = 6. AMP? IN V. 


| 0,80476 

b 0.O1918 

7 0,< 97 40 

11 0.00105 

13 0.00068 

17 0.00038 

19 0.O0025 

FUR RPM s 500.; I DC 
| 3.41081 

5 0.61701 

7 0.00819 

11 0. 00098 

13 O.00072 

17 0,00038 

19 0.60O25 


73.32435 
-31.11028 
-84.51469 
134.53249 
42.42627 
261.08836 
191.31723 

= 6. AMP; IN V 
98.77277 
-31.27879 
-84.35751 
134.71308 
42.24791 
261.73957 
191.68243 


FREQ. -20. HZ 6 WC= 300, 
0.80476 73.32652 

g.pmv 


. .00740 
0.00106 
0.06068 
0.O0O38 
O.O0025 


-31.10732 
-84.54451 
134.53623 
42.55904 
260.81069 
191 .1 3576 


FREQ. =20. HZ & »»C~ 300 


sriiosi 

0.01698 

0*00820 

0.00100 

0.00071 

0.00037 

0.00026 


98.77671 
-31.24336 
-84.44720 
134.49919 
n 42, 58531 
260.79729 
191.14775 




TABHr 3. t> 
* * * * * 


FOR KP« 

1 

s 

7 

11 

13 

17 

19 

Fur rpm 

l 

s 

7 

It 

1.3 

17 

19 

FOR RpM 

1 

5 

7 

11 

13 

17 

19 

FOR RPM : 
1 
5 

7 

11 

13 

17 

19 

FOR RPM 

5 

7 

11 

13 

17 

19 

FUR RPM 
1 

' 5 
7 

11 

it 

19 


VIA FRfiO. nr>MA£f 

>n-y.U I Ooa bubal MAGNITUDE ' PHASE 


= IOC,; toC 

<J , H < 1 4 7 to 
0 , h 1 9 .1 d 
0.< 07 40 
0,00105 
y » t o C 6 8 
0 , Out 3b 
0.0uu?5 


- 6. AMP; Iwv, 
73.3748b 
-31 .1 lo?& 
”84,51 469 
134,53249 
, 42,42627 
261.08836 
191,31723 


F 


= £ >Oto.; me 
3,41 OR 1 
0 ,01701 
6.00419 
0 .00098 
0,09072 
O.inn 3d 

6,00025 

= 500.? IOC 
1 • 3 lo58 
6.01983 
6 ,011 0 8 
0 ,00289 
U.OU182 
6 , 0 6 U 5 to 
9.00035 

= 200 ? IOC 
1,06012 
0.02199 
0.01286 
0,00 4 47 
0.00335 
0 .0017(1 
0.00137 

=. 500„? IDC = 
3,45 333 
6,62313 
6.0158? 

0. 00524 
0.00433 
0 . o 0 2 2 4 
0,00199 


= 500 
. 3 * 4 5 f „ v 
0.07343 
0.01588 
0.00524 
0,00434 
0.00225 
0.00199 


IDC 


=6. AMP; INC. 

98.77277 
”31.27879 
-84.35/51 
l 34.7130 y 
42,24791 
261.73957 
191 ,6824j 

= B.AMPflNV, 

35. 6445o 
22.83890 
" -4,9047 0 ■ 

”62.91205 
266,12543 
194.26089 
148.3 5403 

=4. AMP; INV. 
09,04495 
64 ,46596 
53.38764 
31.O013B 
20.04493 
” 2.42 762 
“13.5210! 

7 6 -AMP? IN tf, 

1 2 2 , / 1 6 6 0 
90.69155 
90.48001 
99,71415 
89,54817 
89.09412 
88.94177 

?22* f 87l| NV 

91.04879 
90.97852 
90.50157 
90,47295 . 
90,|l304 ' 
90.29401 


0.804 / 6 
0.01917 
0.00740 
• 0.00106 
0,00068 
0.00038 
0.00025 

FRFQ =20 
3.41051 
0.01698 
0.00820 
0.00100 
0.00071 
0.00037 
0.00026 


= 20,117 & WC 


3 0 0 


73,32652 
-31.10732 
-84.64451 
134.53023 
42. 559 o 4 
260.81069 
1 91 . 1 1576 

HZ 4, wC= 300. 
98.7/671 
”31 .24336 
-84.44720 
134,49919 
42.585 31 
2o0. 79729 


191,14775 


FRFQ. =30. HZ 6 WC 


” * - ne * ™ ^ 

1.34655 
0.01983 
0.01107 
0,00289 
0.001B2 
0.00066 
0,00035 

FRFy, =15 
i.OoOio 
0.02200 
0.01285 
0.00448 
0.00335 
0.00170 
0.00137 

FKEy, =20. 
3.45320 
0.02347 
0.01586 
0.00526 
0.00433 
0,00226 
0,00198 


800 


85,64724 
22.87026 
-4.92045 
”«2.9y4to6 
266.05194 
1 93,78426 
148.23648 

.HZ 6 WC= 1000, 
99.05103 
to 4 . 4 8 3 y 3 
53,38785 
3 1 « 0 u 9 8 7 

”2*2?3d? 

”13.53293 

HZ 6 WC=9999 9 . 
122.72068 
90.69465 
90.49350 
89.72814 
89.56972 
89.1212? 
88.970/0 


FREQ, =20 .HZ WC=999999< 
122.79196 


*:o!13V 

6100526 

0.00433 

0,00226 

0.00198 


91.05105 

90.99247 

90.51206 

90,49623 

90.33263 

90.32490 
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CHAPTER 


4 


computation 


OF THE electro; iagketic torque of 

CURREI IT FED INDUCTION MOTOR 


THE 


4 . 1 INTRODUCTION 


The expression for the instantaneous electromagnetic 
torque in terms of dq stator ant rotor currents is given by 
the equation (1.8). This equation, together with the expre- 
ssions of dq stator and rotor currents, has been used to 
compute the torque produced by the squirrel cage induction 
motor fed by a three phase current source inverter. 

In the previous chanters, both the time domain and the 
frequency domain expressions for the dq stator and rotor 
currents were obtained. The time domain expressions obtained 
for the dq stator and rotor currents can be directly used in 
the equation (1.8) to obtain the instantaneous torque waveform 

The torque spectra .can also be obtained using the harmoni< 
component values of the stator and rotor currents in equation , 
(1.3). The instantaneous torque expression is given by the 
terms, which are product of the stator and rotor instantaneous 
currents. Thus if a nth harmonic of the stator current and an' : 
nth harmonic of the rotor current, referred to the stator is ; 
assumed then these will produce (m-i-n)th and the (m-n)th harmoni 
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component of the torcue. It has been shown in this chapter 
that the contribution to the torque is only because of one 
of these terms. Further it has been shown that the torque 
harmonic components are multiples of three for the case of ge 
ral stator current. For the case of motor fed by the curren 
source inverter , the torque harmonic frequencies are multiple 
of six times the inverter frequency. 

A torque harmonic has contributions from the various 
combinations of the stator and rotor current harmonics. But 
as the frequency of the harmonic current increases the contri- 
bution of these to the torque harmonic goes on decreasing. Ii 
facx the contribution due to two dominant harmonics of curren - 
is seen to aporoximate fairly well to the actual value of the 
torque harmonic. 

In Sectio n 4.2 the nature of the harmonic components of 
the torque in any arbitary reference frame is discussed. The 
torque harmonic frequencies are found to multiple of three 
times the fundamental frequency for the case of general stator 
current. The general expression for the torque produced due 
to nth harmonic of stator current and mth harmonic of rotor 
current, when referred to the rotor, has been computed in this 
section. 

Sec tio n 4,3 gives the procedures for the computation of 
the torque. One uses the instantaneous stator and rotor 
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current expressions to obtain the electromagnetic torque 
produced. In the other the harmonic components of the torque 
have been computed using the harmonic component values of the 
stator and rotor current. 

In Sec bion 4.4 the stator current waveform of Fig. 2.^- he 
been used to compute the torque produced, using both the pro- 
cedures of Section 4.3. The harmonic components of the instai 
taneous torque waveform obtained are calculated. These are 
compared with the torque harmonic values obtained directly 
through the other procedure. This verifies that both are 
procedures compute the same result. 

4.2 NATURE OF THE HARi IONIC COiTPONENTS INI THE TORQUE 

In this section, the nature of the harmonic components 
in the torque produced by the squirrel cage induction motor 
has been established. To start with a arbitary stator current 
waveform is studied. The expression for torque harmonic pro- 
duced due to the interaction of nth harmonic of stator current 
and mth harmonic of the rotor current, when referred to the 
stator, is computed, in a arbitary reference frame. 

It turns out that torque harmonic is independent of the 
choice of the reference frame and the torque harmonics have 
frequencies which are multiples of three times the stator 
current frequency. 
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Th; stator current of the motor when fed by the current 
source inverter is symmetric about oat = n. Also it does not 
have tviplens. The nature of the torque for such a case of 
the stator current is also discussed in this section. Torque 
harmonic freauencies are multiples of six times the inverter 
frequency in this case. 

4.2.1 Torque harmonic in an arbitary reference frame 


The torque harmonic produced by the squirrel cage induc- 
tion motor fed by a three phase periodic nonsinusoiclal current 
of frequency u, has. been computed in this section. The 
stator current harmonics will have the frequencies of the 
typo ! nio ! , whore n is a positive integer. The rotor currents 
referred to stator are computed after applying a transfer 
function (Fig. 2.2) to the stator currents. Thus the fre- 
quencies of the rotor currents referred to' the stator are als< 
of the type ! mw'. 


The torque produced due to interaction of a stator 
current of frequency 1 nto 7 and a rotor current of frequency 
1 mco’j when referred to the stator has boon studied in this 
section. The rotor is assumed to be rotating with a constant 
speed in the direction of field produced by fundamental 
stator current. Then the harmonics with value of ! m 1 of the 
type (3p--l) produce the flux, revolving in opposite direction 
to the rotor speed. Hero, p = 1,2, 3, 4 ... 
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Thus the actual rotor currents will have the frequencies of 
the type (moo + oo ) where positive sign is taken for the case 
v/her o 'm 5 is of the type (3p-l) , and negative otherwise. 

The expressions for the stator and rotor currents with a 
arbitarv values of magnitudes and phases, I , a and I_, a 




5* 

5 

X X 

respectively, can b 

>e written 

as 



^al ~ X s 

cos[no)t - 

a } 

s J 


(4.1) 

i bl X s 

cos[n(’oot ■ 

„ ?-E\ 1 

3 ' cu s J 


(4.2) 

hi =I s 

cosLn(cot ■ 

4% \ ^ 

“ 3") “ K s j 


(4.3) 

ha = z x 

COS [moot ;r 

“r* “ a r ] 


(4.4) 

^b 2 ~ ^r 

cos[m(oot • 

“ T 1 ) + “V fc “ a r l 


(4.5) 

i c 2 = I r 

cos [m( tot • 

“ 3 ^) ± w^t - a r ] 


(4.6) 


where lower sign is for the case when m is of the type 
(3p-l). Upper sign is fox' all other cases. 

The currents in an arbitary reference frame [3] are related to 
lino currents i a 2 ^ ^ ^tc., from equations (l.l) to 

( 1.4) , as 
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idl = ^ 2 / 3 ) 

i oI = (2/3) 
i e2 = (2/3) 

V = ( ? -/3) 


[i a i cos© + i> ol cos(9 - ) + cos(© + §“)] ( A.l 
[“i al sin© “ i^ sin(©-*|~) “ i ^ sin(© + §“)] (4.8 
[i a 2 C0S I 3 + cos(p - |p-) + i c2 cos(p + |^)] (4.9 
[”i a2 sin3 - i-^ sin(,3 - §^)"i c 2 sin( 3 + pp) ] (4.10 


,3 = © - © (4.11 

r 

©^ = u t (4.12 


In above w is the constant rotor speed. (o©/dt) is the speed 
of the arbitary reference frame. 

The torque is related, to clq currents from equation (1.3) as 


T — 


ii [ i , 
c L ql 


•d2 


■dl 


V 


(4.14; 


where 


H c = M( m/2) (P/2) , . (4.15) 

here m is the number of phase and P the number of poles. 

Substituting for dq currents in equation (4.14) from 
equation (4.7) to oqn. (4.10) and simplifying gives 
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h 

= (*V9) 

[ c [e l ( 

t) 

sin 

(p- 

G) 

+ e, 

2 (t) 

s: 

Ln(3- 

-9- §“) 

+ 






+ e 3 ( 

t) 

sin 

(p 

-0+ 

2tx 
' 3 

)] 





(4, 

,16) 

ah 2 

i re 

















X 

(t) = 

(i 

al 

X a2 

T- 

1 bl 

1 b2 

~h 

1 cl 

i c2 ) 


(4, 

,17) 



e 2 

(t) = 

(i 

al 

X b2 

*r 

1 bl 

X c2 

4* 

1 cl 

i a2 ) 


(4, 

.18) 



e Q 

(t) = 

(i 

a2 

1 bl 

X. 

i b2 

1 cl 

+ 

i c2 

i al^ > 


(4, 

.19) 

Sub 

S titu C3 

.ng 

for i 

al^ 

t) , 

ia: 

0 ( t 

) in equal 

;ion 

(4.17) 

f rom 



equ 

.? tions 


1) to 

(4 

.6) 

. 0 

ive 

s 









e^(t) = I g I r [ cos( nwt“Oc s ) cos(mcot + - a ) + 

+ cos(a(cot*“ §~) -a ) cos(m(wt-~ ?r-) + co t - a ) + 

v-> b O ^ J X X 

+ cos(n(wt- *p-) -c: s ) cos(ui(wt“ ^~) + io r t-c: )] 

This implies 

e j_(t) = (l s I r /2)[cos [(n+m)a)t~(a s +a^)+oj r t] + 

r\ a 

+ cos [ (n+m)(oJt- •j^)+ a) r t~(a Q +a r )] +cos[ (n+m) (u)t~ 3 ™) 

± w r t~(a s +a r )] -r cos [(n-m)wt + u r t“(cs s -a r )] + 

+ cos f(n-m)(tot- v~) - w_t~(o: -a)! + 

+ cos [ (n-m) (u>t~ 


( 4 . 20 ) 
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ihis can be shown to simplif y to 


Ojit) = (3I s I ;: /z) [p s cos[ (n+n)ut+w r t-(a s +a r ) ] + 
J " P C 1 cos [ (n-m)wt + a) r t - ( C '- S ~ C ' T ) ] ] 


(4 .21) 


P< 


if (n-<-m) is of type ’ 3p’ 


0 


P' = i if (n-m) is of type ' 3p ! 


>therwise 


(4.22) 


(4.23) 


Here p is any integer. 

Similarly from equations (4.10) and (4.19) s it can be shown 


C i~13 


&o(t) = (31 I r /2) [p cos [(n+m)cot - + to t -(a +a )] + 


'2 


s r J 


+ p c , cos [(n-m)tot + + to_.pt - ( a s “ a r )] 1 (4.24) 

,( t) - (3I s I r /2) [p s cos [(n-Hn)cot+ ^ 


2Tcm + p , 

" l,”( a s +a r ) J + 


r p., cos [( n-m) cot - -5— + ta_,t-(c: ~o: )] ] (4.25) 

t-’ O JL ip X* 


The expressions of e^t) , © 2 (t) anr ' 1 ^(t) t row. equations 
(4.21) s (4.24) end (4.25) are substituted in equation (4.16). 
Table 4.1 lists the final expression of torque obtained after 
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simplification for various combinations of ! n' and ! m' values. 
In the Table 4.1, ®n* corresponds to stator current harmonic 
and ’ m® the rotor current harmonic when referred to the stator 
and p is a positive integer, i.e. 0,1,2, 3 ... 

It can he seen from Table 4.1 that the torque spectrum 
has frequencies which are multiple of ’ 3m ’ . 

The frequency components of the type (3p+2) produce the 
flux rotating In opposite direction to the flux produced by the 
fundamental. The different torque expressions written in 
Table 4.1 can be written as one general expression if the 
frequency component of type (3p+2) and its phase is assigned 
a negative value. This implies that m or n takes the values as 


-2,3,4,- 


(4.2b ) 


If this is done, the torque is given as 

T = 21 I M sin[(n“-m)wt“(a; , '-a' ) ] if (n-m) is multiple of -3 

CJ S _l C Si 1 

= 0 otherwise (4. on) 

where n and m takes the values as given by equation (4.27) and 


*>c 

a s ~ a s if n is positive 

= if n i$ negative (4. 4*0 

c:^ = a r if m is positive 


if m is negative 


(4.TSU 
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Thus if n and m are represented by values given by equation 

(4.27) s it is seen that the torque harmonic given by the 
difference of the stator and rotor current frequencies is 
produced. The equation for the torque (4.28) can be seen to 
be independent of the choice of the reference frame. 

4.2.2 Harmonic components of the torque for the induction 
motor fed by a current source inverter 

The stator current of the induction motor fed by a current 
source inverter is antisymmetric about tot = it;, i.e. i(9) = 

-i(0 + it). This implies the absence of the even harmonics 
in the stator current. The triplens are also absent because 
the sum of three phase current will be zero in this case. 

Thus the permissible values of stator harmonics from equation 

(4.27) are 

1, -5, 7, -11,13, “17, 19 (4.30) 

In a induction motor rotating at constant speed the induced 
rotor currents, referred to the stator will also have harmonics 
given by equation (4.30). Thus stator and rotor currents have 
odd and non triplens values only. 

It has been shown in Sec. 4.2.1 that torquo has contribu- 
tion from the harmonics when (n-m) is triplen. Values of 1 n* 
and ’m’ in this case, equation (4.30), are odd, and thus 
(n-m) is even. Thus the torque spectrum of the induction motor 
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using fed by the current source inverter have harmonics which 
arc multiplies of six. 

In previous chapters, the harmonics of the dq frame stator 
and the rotor currents have been obtained. The nth harmonic 
an:’ nth harmonic of stator and rotor d axis current can be 


i C !l = cos ( n0j t - a^) (4.31) 

i ^2 = £9 cos(mot - cc 0 ) ' (4.32) 


I 1 ,I^ ) ,c. 1 and cm can be related to 1 , t , c: and as follows. 
Substi taring for the values for the values of current in equa-- 
tion (4.' 7 ) from equations (4.1) to (4.3). 

rs 

iji = -v I^[cos(ni)t-r^) cos© + 

C J. O o !o 

[cos [n(iot - §~)~c , s ) ] cos(© - *|“)] + 

-i- [cos[n(mt - &)-g ) ] cos(© - P ) ] 

u O s J O 

Simplifying this, v:e obtain 

i^-i = I s cos(no)t - c: s + 9) (4.33) 

’/■hero upper sign is v;hen n is of type (6p~l) and lower 
o ';hs r’ A s e . 

Comparing equations (4.33) and (4.31) gives 
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X s = I 1 


a. = a, +9 
s 1 “ 


Similarly it can sbe shown that 


t = t 
i r J "2 


( & 2 i P i w^t) 


= CCq + 0 


(4.34) 

(4 .35) 

( 4.36) 

(4.37) 
(4.37b) 


>u 


bstitating fron equations (4.34) to (4.37) in equation (4. 28) 


= 21-^2 i.l c sin[ (n-m)a)t - (c^-o^)] (4.38) 

Here n and m takes value as given in equation (4.31). 

The signs of a ^ anc] a 2 are reversed if n or u respectively has 
a negative value. 


It should be further noted that the harmonics of the 
pseudo rotor currents, X^ and X 2 * rather than rotor currents 
have been computed in the previous chapters. The dependence 
of torque on rotor currents can be changed to pseudo rotor 
currents by equations (2.42) and (2.43). Substituting these 
in equation (4.14) gives 


T = M, 


^V X 2 


22 


-dl' 


■dl 


X, 


M 

L 


22 


"ql 


)] 
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This implies that 


T 


“ l c L ^ql 


An 



(4.39) 


It is seen that this equation (4.39) is similar to the equa- 
tion (4.1) except for the change of i^ ancl iqp ' i:o ^2 anc * ^1* 
respectively. and X-, can be seen as ’ d ! and q axis compo- 

nent from equations (2.42) and (2.43). Thus a equation similar 
to equation (4.33) can be obtained in terms of the harmonic 
components of the escudo rotor currents as 


'T* 

i 


q 


; a 1 is 
values 


= 21 ) X x ; ,; c sin [(n-:.)wt-(c 1 -c: x )] 

is mth harmonic magnitude of X 0 . a x is 
positive Slid, negative of phase otherwise 
as given in equation (4.31). 

THODS FOR THE COMPUTATION OF TIE TORpUE 


(4.40) 


its phase if 
’ m* take 


The expression relating the dq stator and rotor currents 
to the instantaneous electromagnetic torque is given by equa- 
tion (*<5.i). In previous chapters, two methods for the computa- 
tion of the rotor currents from a given stator current have been 
described. 

The first method results in the time domain solution of the 
rotor current and the other gives the harmonic components of the 
rotor current. 
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A method to compute the torque produced by the squirrel 
cage induction is to substitute this time domain solution of 
the rotor current obtained, in equation (4.1) to get the time 
domain torque expression. This method is described in 
Sec . 4 .3.1. 

Other method uses the values of the harmonic components of 
stator and rotor current obtained. Equation (4.39) gives the 
expression for the harmonic torque produced due to the interac- 
tion of the nth harmonic of the stator and mth harmonic of the 
rotor current. For the system of motor being fed by the 
current source inverter, m or n takes values as given by the 
equation (4.31). For comoutation of particular torque harmonic 
component, which is multiple of six (Sec. 4.2.2), all the 
possible combinations of the permissible values of * m* and 'n' 
are taken which contribute to this torque harmonic. Their 
contributions to the torque can be obtained from equation 
(4.38). In Section 4.3.2, the procedure to evaluate all the 
possible combinations for a particular harmonic has been men- 
tioned 4 

4.3.1 Time domain solution of the torque 

In previous chapters, the method for the time domain 
solution of rotor currents has been given. There the cases of 
stationary and rotating rotors have been dealt with separately. 
For the case of stationary rotor, the rotor currents expre- 
ssions have been obtained. This can be substituted in 
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equation (4.14) to obtain the torque produced for the stationary 
rotor case. 

For the case of rotating rotor the solutions obtained for 
and X 2 can be substituted in equation (4.40) to obtain the 
ins tantaneous torque values. 

It has been shown in Section 4.2.2 that the torque has 
harmonics to be multiple of six times the inverter frequency. 

It implies that the torque waveform repeats after every 60° 
interval. Thus the computation for this section has to be 
done only for any one of six intervals of the inverter. 

4 - 1-2 frequency '^io ' So\utio o of "fKe. "Wityxe. 

The value of a particular torque harmonic can be computed 
through equation (4.40). For this computation, first of all, 
the possible combination of tho current harmonics of the stator 
and rotor have to be obtained which contribute to this torque 
harmonic . 

The stator and rotor currents have infinite number of the 
harmonic components. The magnitudes of these harmonic go on 
decreasing as the number of harmonic increases. Thus though 
there are infinite number of the combination which will give 

the particular torque harmonic, the contribution because of 

1 

higher order current harmonics will be very insignificant. 

In the section the evaluations of stator and pseudo rotor 
currents harmonic combinations is done in a special sequence. 
This sequence has a property that at each step all the 
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permissible combination of harmonic currents up to a specific 
frequency are computed. The torque value is updated with 
every new combination taken till it is observed that the 
contribution because of coming new pairs to the torque is 
less than a limit (say 0.001^) of the last updated value. 

This calculation of the torque harmonic has the following 
four steps. 

(i) To got the initial pair of stator and rotor frequency to 
start the torque computation. This has been shown in Fig. 4.1. 

(ii) To compute the torque value by a given pair of frequencies. 
This computation has been shown in Fig. 4.2. 

(iii) To check when to terminate the present process of compu- 
tation. This is done by checking the contribution of the last 
pair to the torque; Further calculations end if this contri- 
bution is less than a limit (say O.OOl/O and the two possible 
combinations with a component of pair 1 had been considered. 

If this is not the case we move to next step. 

(iv) This step computes the next permissible pair which con- 
tributes to the chosen torque harmonic. 

The flow chart of the procedure for computation of the 
chosen torque harmonic has been given in Fig. 4.3. 
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4.4 CALCULATIONS FOR A KNOWN MACHINE PARAMETERS 

A computer program is developed to compute torque harmonic 
using the flow chart of .Fig* 4,3. The time domain torque Is 
obtained via Sec* 4.3.1 approach. The- listing or? -computer pro- 
gram is given in the Appendix B. 

The average and sixth harmonic torque of induction motor 
produced by the stator current of Fig. 2.1 has been plotted in 
Fig, 4.4 as the function of rotor speed. It is seen that 
sixth harmonic torque is maximum at synchronous speed. 

Fig. 4*5 gives the relative contents of various torque 
harmonics at three rotor speed. It can be noted from this 
f igure that as the rotor speed is increased the torque harmonic 
components become mare and more dominant. 

Fig. 4.6 plots the time domain torque waveform* This wave- 
form has been drawn for 60° interval only because the torque 
waveform repeats after every 60° interval (Sec* 4*2). 

Table 4.2 compares the values of torque harmonics obtained 
from following methods j 

(a) Using Sec. 4,3 method of obtaining torque harmonics from 
harmonics of currents. 

(b) Obtaining the harmonic components of the time domain 
solution of the torque* obtained through Sec. 4.2 approach. 
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(c) Obtaining tho value of torque harmonic by considering 
the dominant current harmonics only. For tho torquo harmonic 
frequency of 6 nw, tho dominant current harmonics have fre- 
quency (6n-l)w and (6n+l)to because those interact with the 
fundamental to produce the harmonic torque. 

It is seen from Table 4.1 that all the three methods give 


tho similar results. 
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START 


Divide the nth torque harmonic component 
value into two halves, n/2 and n/2 




i Go on incrementing one value and decre- ! 
I menting the other till neither i s even or' 
j triplen ! 
I | 


The number corresponding to type (6p-l) 
is assigned the negative value 



The two numbers form the initial pair of ' 
frequencies 1 j 
! 


V 


END ; 


Fig. 4.1 s Step_l Com putation of jt he in itda' 1 

s tar ting jp air rogue n cyf^f-~th^s t a to r 

and rotor current — 
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start' 


Assign a boolean variable (say DONEKEV)] 
a value true ! 




j Suopcse the input pair is (n,m) computes 
; nth harmonic of stator and mth harmonic of 
pseudo rotor current 


J Compute torque value using equation (4.41) | 


| Add this computed torque value to the previous | 
j value keeping care of phase as this is vector | 
i addition I 

i i 



Toggle variable DONEREV 


i 

! 


Reverse the 
pair to 
(m,n) 


Fiu'.o 4.2 ; Step 2 U pdating of torque with the 

new pa ir v alue 


No 


/ 


X 


- 'MS DON ERE V> 


True?/ 

i 




f 


Yes 


END . 
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( START : 


(i) 


A boolean variable DONEPLUS 
is assigned true if chosen 
torque frequency divided by 
six is odd else is assigned 
pulse 


(ii) DIFVAR 


- 6 


| Step (i) of 
; Fig . 4.1 


| Increment 1 
1 DIFVAR by S 
| number 6 j 

‘ 7 — IT 

! j 


Increment | ! 


Decrement 

pair numbers j 1 • 

< 

pair numbers 

by DIFVAR j 




by DIFVAR 

-r-t 
* ( ' 


Step (ii)' 
[of- Fig. 4. 2 




No 



IS \ 
DONE PLUS " 
true? 


Yes 


\ 




/step (iii) >- 

/ 

/ 


Toggle 
DONE PLUS 


\ 


. continue 


Terminate 


\ END 


Fig, 4 . 3 Fl ow cha rt for the computation of a 
pa rticular to rqu e harm onic 
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CHAPTER 5 

REDUCTION OF TORQUE HARMONICS BY MODULATION OF D.C. 

INPUT CURRENT 


5 . 1 INTRODUCTION 

The induction motor being fed by a current source 
inverter produces harmonic torques. The harmonic torques 
can be reduced by modification of the stator current waveform 
This can be done either by modulating the dc current input of 
the inverter or by modulation within the inverter. In this 
chapter the former technique of modulation has been studied. 

The latter technique is referred as pulse width modulation in 
the literature [10], 

The aim of the modulation is to reduce the harmonic 
torques for the same average torque. The interest Is directed 
more toward reducing the dominant torque harmonics, 6th and 
12th. The higher order harmonics apart from being lower in 
magnitude are not reflected in the shaft speed because of the : 
large inertia of the mechanical system. 

The reduction in the torque harmonic can be noted by 

the increase in factor (T aV /T n ) after modulation, where 1 T av ’ ; 

and 'T ' refer to average and nth harmonic torques respecti- 

n t i 

vely. This is so because torque can be seen to be proportional, 
to the square of stator current, from equation (4.33). Thus 
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the inverter input current when multiplied by the factor 
[(T n )no mod/( T n )mod] gives the average torque after 
modulation equal to case without modulation. This allows 
us to compare the relative contents of harmonic torques 
ox two torque spectrum of different average torques. 

It has been shown that the modulating waveform has to 
fulfill a general requirement so as to produce a balanced 
three phase stator current. This modulating current should 
repeat after every 60° interval, that is, the modulating 
current Jrrequency should be six times the inverter frequency. 

In this chapter the cases of modulating current waveform 
of exponential or of consinusoidal nature have been studied. 
Having assumed a modulating current, satisfying the above 
requirements, the torque produced by the motor has been com- 
puted both in time domain and frequency domain (Sec. 4.3). 

Time domain solution is obtained by first evaluating 
pseudo rotor currents using method in Section 2.2. These 
expression when substituted in equation (4.39) gives the 
time domain torque expression. 

For computation of the torque spectrum in frequency 
domain, the harmonic components of the stator dq current is 
obtained through a computer program. The harmonic components 
of pseudo rotor currents are evaluated by proceeding as in 
Section 2.3. Then procedure in Sec. 4.3.2 can be used to 
compute the torque spectrum. 
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In Section 5.2 the origin of the general requirement on 
the modulating waveform has been shown. In Section 5.3 the 
case of exponential waveform as the modulating waveform has 
been studied. In this case the amplitude and the time constant 
of the exponential waveforms are the parameters for the study. 
In Section 5.4 the case of consinusoidal modulating waveform 
has been studied. In this case, the magnitude, frequency and 
the phase of the modulating waveform act as parameters for 
study. 

It should be noted here that in this chapter, for the 
purpose of simplification of analysis, the current source 
inverter has been assumed to be ideal. 

5.2 GENERAL REQUIREMENTS OF MODULATING WAVEFORM 

The current source inverter produces balanced three phase 
line currents. These currents have following property : 

(i) They are displaced from one another by 120° 

(ii) Each of the line current is antisymmetric (i.e., 

i(©) = -i(it + ©)) . 

Let the current during interval I in phase ’a’, be a general 
function of time f(t), where t refers to the time with origin 
at the start of an interval under consideration. For example, 
f(t) may be a exponential waveform as shown in Fig. 5.1. Then 
the intervals III and V of phases 1 b' and 1 c* respectively will 
also have the same current f(t ) due to property (i) above. 
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Interval IV phase ’a' current is - f(t) due to property 
(ii) above. Thus phase 'b* and 1 c f have currents of type 
- f(t) during intervals VI and II respectively. 

Thus during all six intervals, the conducting phases 
carry the same function of current f(t). Therefore, inver- 
ter input current should be same during all intervals to 
satisfy the above condition. This implies that the modulating 
current should repeat after every 60° interval, that is, 
it should have a frequency equal to six times the inverter 
frequency. 

5.3 EXPONENTIAL MODULATION 

The general expression for the exponentially modulated 
inverter input d.c current, i(t), during a 60° interval of 
inverter can be written as 

i(t) = I + K(e^ t - 1) (5.1) 

whore ! K ! and ! j3 ' are the magnitude and inverse of time con- 
stant of the modulating waveform. * I’ is the unmodulated input 
d.c. current. Here the origin of time, t, is at the start of 
each of the intervals. The three phase currents obtained by 
this modulated inverter input d.c. current are drawn in , 

Fig. 5.1. To study the effect of the modulation on torque 
harmonics, both time and frequency domain solutions for the 
torque are obtained. 



From the time domain solution of the torque it has 
been shown in this section that the torque harmonics can 
be reduced to zero for the case of stationary rotor, if 
K and (3 are chosen to be * I* and * a' (eqn, (2,18)) res- 
pectively. For the rotating rotor case it has been proven 
that it is not possible to get a constant average torque by 
using the exponential modulation. 

The effects of variations in K and p on the harmonic 
torque spectrum have been studied through frequency domain 
torque computation technique. The results obtain show 
that any exponential modulation improves the performance index, 
T av /T n s for all rotor speeds below the synchronous speed. But 
in case of speeds greater then synchronous speed, there is a 
deterioration in this performance index if due to any exponen- 
tial modulation. 

5.3.1 Time domain analysis for the exponential modulation 

In this section, the analysis for the stationary and 
rotating rotor cases has been done separately. It has been 
shown that a constant average torque can be obtained for the 
case of stationary rotor but not dyring rotating rotor by 
exponential modulation. 

5.3. 1.1 Case of stationary rotor 

During interval I, Fig, 5.1, the currents are 
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where 

These 


i = K e 
a 


pt + j* 


i h = -(K e Pt + I*) 


i c =° 




I = (I - K) 


can be transformed to dq frame as 

i di = [K e pt + I*] 

V = - f [K oP* + I*] 
qX f3 


(5.2) 

(5.3) 

(5.4) 

(5.5) 

(5.6) 

(5.7) 


The time domain solution for the torque during this interval I 
can be obtained using equation (4.40). For this X^(t) and 
X^t) have to bo calculated. 

For the case of stationary rotor, i.e., oo r = 0, equation 
(2.51) gives 

pX 2 = -aX 2 + ai ql 
This implies 


(a+p) Xg - ai d 

Substituting for i^^ i n this equation from equation (5.6), 
gi. ves 

(a+p) X 2 = y ~ — ■* a[K e^ + I ] (5.8) 

z L 22 
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Tho solution of this equation (5.8) is 
X 2 (t) - C 1 e~ at + A ] _G^ t + B 1 

where 

1C, K 
A 1 ” a+p~ 

* 


(5.9) 


(5.10) 

(5.11) 


K 9 =f~a ' (5.12) 

z l 22 

and Cj| is a constant. 

Similarly, from equation (2.50) we obtain 

Xf(t) = C 2 e~ at + D x + E 1 (5.13) 

where 

1 *2 K 
D i - - p • T5+f) 

S i* 

b** 32; mm «3*=««*a* fa trsttux* 

1 f3 a 

and C 2 is a constant. 

Substituting for i^, ^qi» ^1 an< ^ ^2 ^ rom a 9 ua tions (5*6), 
(5,7), (5.14) and (5.9) into equation (4.40), we obtain tho 
torque expression during interval I, after simplifying as 

~ at + K e (P ~ a)t ] 


(5.14) 

(5.15) 


T q = " M c^ C 2 + T§ } Q 
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The time dependence of torque in above equation can be removed 
if K and p are chosen as 

K = I , p = a (5.16) 


For this choice* 


T 


q 





I 


(5.17) 


that is * a constant torque is obtained. 


The values of and Cl, depend upon the initial conditions 
and can be computed by the analysis, as shown in Sec. 2.2. 
Using the method of Section 2,2, the following values can be 
obtained 


'1 1 


c 2 | 


V3/2 


( 4 ^+ 3 ) | 


! -1T3/2 
L 




L Y*3 


i ^ X 2 +X P 


x* = (£ - e “* a/3a) ) 


(5.18) 


* 





+ B 


1 


x 3 “ ^1 + ®1 


(5.19) 
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Table 5.1 lists the expressions of i^, i ^ and 

T with values of K and p from equation (5.16) for intervals 

v-3 

II and III. Here C^, C^, and are constants which depend 
on the initial conditions and can be calculated as in 

Section 2.2. In this, t 2 and tg are as defined in Section 2.2. 

The torque is seen to be constant for intervals II and III for 

K,p choice of equation (5.16). Intervals IV, V and VI are 

similar to intervals I, II and III respectively because of 
symmetry of current about cot = ir« 

Substituting, from equation (5.16), in equation (5.1), 

gives 

i(t) = I e at (5.20) 

This equation (5.20) gives the form of the input current 
during an interval, which will produce a constant torque for 
the stationary rotor case. 

5. 3. 1.2 Case of rotating rotor 

The pseudo rotor currents are related to stator da 
currents from equations (2.52) and (2.53) as 

(p 2 +2ap+a 2 ^)X 1 = K ia l ql + K lWr i dl + K lP i ql (5.21) 

(p +2ap+a +w r ) X 2 — ~ (5*22) 

Let us confute the torque expression for a period during inter- 
val I i.e. for f u)t‘ values between 0 and ir/3. 
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Using values of and i ^ for this period in (5,21) 
and (5.22) from equations (5.6) and (5.7), the expression 
for and y^_ are are obtained as 


X 2 
where 


= C ^ e cos(w r t + 0 ^) 

B "t 

4- A^2 & 4* 

(5.23) 

= Cg e~ at cos(m r t + 0 2 ) 

+ A 3 + Bg 

(5.24) 


K q x 

^ +2ap+a^+m'‘ 


(5.25) 


q x = K 1 (w r - ^--1) 


B 


q 9 


* 


kjssskxs. jgr - n . 

2 (a 2 -K»p 


q 2 = Kl^r - p) 



( p +2ap+a +wp 


0 ) 


q 3 = + a + p) 


B 


q 4 


* 


3 

°*4 


(T^f> 
K l^f§ + a > 


(5.26) 

(5.27) 

(5.28) 

(5.29) 

(5.30) 

(5.31) 

(5.32) 
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C 7’ C 8’01 and ^2 in equations (5.23) and (5.24) are constns 
which depend on the initial conditions and can be computed 
by the analysis as done in Section 2.2. 

The expressions for i dl , i ql , Xj _ and Jt, are substituted 
from equations (5.6), (5.7), (5.23) and (5.24) into equation 

(4.40) to get the expression for torque. After siroplif ication 
this equation is given as 

Tq = “ M C ^ K e ^ P c 9 cos(a> r t+0 3 ) + KA 4 e 2 ^ + 

+ (KB 4 +A 4 I )e^ t +C 9 I e at cos(w r t+0 3 )+B 4 I*] (5.33) 


v/here 


A 4 - $ 


B 4 = 5" 


4 K l U) r I 




(a^7, 


and Cqi 03 are obtained from relation 


(5.34) 


(5.35) 


C 

Cg cos(co r t+0 3 ) = C rj cos(w r t+0q) + cos(m r t+0 2 ) (5.36) 

That is f Cg and 0 3 depend upon the initial conditions. : 

To remove the time dependence of equation (5.33) we have the 

choice of two parameters K and p. Sinco there is a term of 
23 1 

typo ' G ’ > to make equation (5.33) time independent, it’s 
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coefficient should be made zero by choice of K or p, This 
is possible only either K or is zero. Since we are 
oealing with rotating rotor case, ^ 0, For exponential 
modulation K ^,0. Thus (K.A 4 ), the coefficient of e 2pt cannot 
be reduced to zero. 

Thus it is not possible to get an exactly constant torque 
for the rotating rotor case, using the exponentially modulated 
waveform. 

Though it is not possible to make equation (5.33) inde~ 
pendent of time with K and j3 as parameters, but it is evident 
from this equation that certain choices of K and p may give 
lower time variations in torque, then others. This can be 
seen in a better way through torque spectrum in frequency 
domain. This has been studied in detail in next section 
through the frequency domain analysis of the torque spectrum 
duo to the current waveform of Figure 5.1. 

5.3.2 Frequency domain analysis for the exponential modulation 

The effects on rhe torque spectrum due to the variations 
in the exponentially modulated waveform, Fig. 5.1, have been 
studied in this section. The torque spectrum produced by this 
stator current waveform has been computed using the method of 
Sec . 4.3.2. 

It is seen from equation (5,1), that for exponential 
modulating waveform there are two control parameters, p and K. 
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In this section the effects on torque spectrum due to the 
variation in the values of p and K of the modulating waveform 
has been studied for various values of rotor speeds. When the 
effects due to variations in ’ p* are studied, ’K* is assigned 
a value equal to ’I*, which gives no torque harmonics for 
stationary rotor case, (Sec. 5. 3. 1.1). Similarly, when effects 
on torque because of variations in ’ K* is studied, 1 p* is 
assigned the value * a f , as this value of ' P ' removes torque 
harmonics for stationary rotor case. Sec. 5. 3. 1*1. 

5. 3. 2.1 Analysis with the variations in the time constant of 
the modulating exponential waveform 

To study the effect on the torque spectrum due to variation 
of P, the inverse of modulating waveform time constant, the 
amplitude of this waveform, K, is taken as I. It is noted that 
this value of K and p equal to * a 1 give a constant torque for 
stationary rotor case. For this choice of K, the expression for 
the modulated inverter input current is given by, from equation 
(5.1) as 

i(t)=Ie ?t (5.21) 

For the current the torque spectrumccan be obtained for 
various values of p , at various rotor speeds. Suppose we are 
interested in reduction of the 6th harmonic of the torque. 

Fig. 5.2 plots the performance index (T av /T^) as function of the 
rotor speed. refers to 6th harmonic torque and T ay to the 
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average torque. In Fig. 5.2, (3, the inverse of the time 
constant of the modulating waveform has been taken as para- 
meter. The case (3 = o, in Fig. 5.2 corresponds to the no 
modulation case. It is evident from this figure that 

(i) For the rotor speeds less then the synchronous speed 

of the motor there is improvement in the performance with aaf* 
value of ' p* of the modulating waveform as compared to the no 
modulation case. As against this, for the rotor speeds 
greater than the synchronous speed of the motor there is no 
improvement with the exponential modulation, rather there is ar~" 
deterioration in the performance. 

(ii) It is seen from Fig. 5.2 that there is a optimum .value of 
P, at particular value of rotor speed which gives the best 
performance. For rotor speeds loss then synchronous speeds, 
this optimum value of {3 increases as the rotor speed increases. 

(iii) (a) The performance index (T^/T^) is a sharp function 

of {3 at rotor- speeds much less, then synchronous speed. However, 
at higher speeds, that is, closer to synchronous speed (T qv /T^ 
is not very sensitive bo values of p. For example ; at rotor 
speed of 300 rpm, ( T av / T 6^ ^ or P ” a 20.3 and for (p=2a) is 
31.4. Similarly at 200 rpm when p changes from a to O , 
(Tqv/T 6) rQC * uces 36.0 to 17.4. At higher speed say at 

560 rpm, (T^/T^) for P = 12a and 16a is 5.87 and 5.06 res- 
pectively. At 580 rpm, (T ov /T^) changes from 1,41 to 1.98 
for change of p from 12a to 16a. 
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( b) In view of above, it is possible to keep p constant in 
the normal operation region of induction motor, i.e. between 
approximately 90/. to 100^ of the synchronous speed and still 
obtain near optimum value of (T /T^) ♦ 

(iv) As the rotor speed tends to zero the optimum value of 
’p 1 approaches the value ’a', as obtained in Sec. 5.3.1. 

5. 3. 2. 2 Analysis with the variations in the amplitude of 
the modulating exponential -waveform 

To study the effects on torque spectrum due to tho varia- 
tions in the amplitude of exponential modulation, K, the value 
of p is kept constant. This constant value is taken as ’a*, 
as this choice of p gives constant torque for stationary rotor 
(Sec. 5. 3. 1.1). In this section, the interest has boon focussec 
on the reduction of the sixth harmonic torque with variation 
in K. 

Fig. 5.3 plots the performance index (T /Tg) as the 
function of rotor speed, varying from zero to synchronous 
speeds. The case K = 0 in the figure corresponds to the case 
with no modulation. It is evident from this figure that the 
effects on torque spectrum due to variations in K are similar 
to • the variations in torque spectrum with p as a parameter. 

That is, for every value of rotor speed, there is an optimum 
choice of the value of K, which gives us the maximum perfor- 
mance index, (T / T^). Similarly, this index is a strong 

3 v v 
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function of ’ K* at rotor speeds much less then synchronous 
speeds. However at higher speeds, that is, speeds closes to 
the synchronous speeds (T /t a ) is not a very sensitive to 
the values of K. In view of this, it is possible to keep K 
constant in the normal operation region of the induction motor 
and still obtain a noar optimum value of {T av /T^) . 

It can be seen' from Fig. 5.2 that as the rotor speed 
increases, the optimum value of K for the best performance 
increases. For rotor speed approaching zero, this optimum 
value of K approaches the value I, inverter input dc current wit 
no modulation. This is in fact the optimum value of K, which 
is obtained for tho case of stationary rotor (Sec. 5. 3. 1.1). 

5. 3. 2. 3 Optimum pair of (K,p) for a given speed 

It has been observed in Sec. 5. 3. 2,1 that for a given rotor 
speed and a assumed constant value of K, there is a value of p, 
which gives the maximum performance index T av /l^. Similarly, 
Sec. 5. 3. 2. 2 shows that when p is taken as constant, this opti- 
mum property is observed with respect to the value of K as well. 
It has been shown in this section that there is a optimum 
choice of tho pair (K,p) which gives a maximum performance index 
at a speed. 

To observe this , the maximum T av /T^ value is obtained for 
various choices of j3 by varying K. The plots of (T av /T£) max 
as a function of J3 and (K) optimum as a function of : 3 has been 



5.17 


plotted in Fig, 5,2 for the ease of rotor speed » 500 rpm. 

F-rom those plots it can be observed that, 

(i) There is an optimum value of the pair (K,p) which gives 

tho largest (T av /T 6 ) maJ[ value. For this pair, the (T av /T 6 ) maJt 
is very largo and it can be assumed that a near constant tor- 
que is obtained with this choice of (K,p ) pair. 

(ii) As the value of (3 increases, tho value of K which gives 
the maximum value of (T^/T^) for tho chosen p, goes on 
decreasing. 

Fig. 5.3 plots the optimum (K»p) values as the 
function of speed. From this plot it is seen that as tho rotor 
speed increase tho optimum value of p goc-s on increasing. It • 
is equal to * a’ , at zero rotor speed, 

5.4 COSINUSOIDAL MODULATION 

Tho general expression for the cosinusoidally modu- 
lated inverter input dc current, i(t), during a 60° interval 
of invertor can bo written as 

i(t) = I + Am cos (u)t + a ) (5.38) 

whore A jn? w m , <x ra are the amplitude, frequency and tho phase 
angle of the modulating waveform. ’I* refers to unmodulated 
d.c. current. The effects on tho torque harmonics due to such 
a modulation has been studied in this section with both the 
time and frequency domain torque solutions. 
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From the time domain solution of the torque it has shown 
that it is not possible to get a constant average torque with 
a cosinusoidal modulation. The effects of variations in A m , 
u> m and <x m on the harmonic torque spectrum have been studied 
through the frequency domain torque computation technique. 

5.4.1 Time domain analysis for tho cosinusoidal modulation 

In this section, it has been argued that it is not possiblt 
to obtain a constant torque with cosinusoidal modulation. The 
cases of stationary and rotating rotors are dealt together. 

Let us compute the torque expression for a period of 
interval I i.o. for * u> t * values given by 

o < tot < ^ 

For this period, the three phase currents for the dc input 
current given by equation (5.38), are 


X a = X+A m cos <“m t+a m> 

(5.39) 

X b = - tl+A m “ st “m t+a m>] 

(5.40) 

i =0 
c 

(5.41) 

These phase currents can be transformed 

equations (l.9)and(1.10) as 

to dq currents from 

i dl= I+A m cos( “ro t+c ‘m ) 

(5.42) 

V 3 * ~ [ I+A m cos^t+otj] 

(5.43) 
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Tho expressions for and can be obtained by solving 

tions (5.21) and (5.22) with above values of i , and i « 

ql q2 


gives 


oqua- 

This 


X 1 (t) = C 1Q e" at cos(w r t+0 4 )+D 1 cos (u m t+a m )+D 2 sinCw^t^)^ 

(5.44) 

t ) = C 1X e“ ai: cos(w r t+0 5 )+D 3 cos(w ra t+a m )+D 4 sin(w m t+a n )+E 2 

(5.45) 

whore 



<*2 


q 

D 


2 


D 


1 

3 


24 


= Am [q 2 (a^+w r +w m )-2aqu} m ]/E 3 
= -[4a 2 ^+(a 2 -Ho^^) 2 ] 

= K l (w r “ 

_ K 1 w m 

- fT 

= Am[2aq 2 m m +q(a 2 -to 2 -w 2 )]/E 3 
q 2 I 

— ' ~7y yy 

a“+u r 

= Am[q 4 (a 2 +w 2 -w 2 )~2aq 5 w m ]/E 3 
== Kj_( a+ p) 


(5.46) 

(5.47) 

(5.48) 

(5.49) 

(5.50) 

(5.51) 

(5.52) 

(5.53) 


q 5 K x w m 


(5.54) 
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D 4 = Aiti[2aq 4 0J m +q 5 (a 2 4w^-a>^)]/E 3 

r. _ ^4 ^ 

C 9 — - W - -"Ty 

a +wj 


(5.55) 

(5.56) 


C -,0 and arc tho constants which depend on the initial con- 
ditions. These can bo evaluated by tho analysis as done in 
Sec. 2.2. 

Tho expressions for and obtained are substi- 

tuted in equation (4.40) to get the expression for torque, Afte; 
simplification this equation is given as 


T q * -M c [C 12 I c- cos(co r t+0 6 )+lD 5 cos(u> m t+a)+ID 6 $in(w m t+a) 4 
+ C l 2 Me”‘ at cos(u> r t+0 6 ) cos(w m t+a) +(^)M D 5 cos( 2 w m t+ 2 a) + 
D 6 sin( 2 w m t + 2a) + (^)M (D 5 +D 6 )] ( 5 . 57 ) 


where 

(5.58) 

(5.59) 


C-j _ 0 and 0 £ are obtained from equation 

C 12 cos(o) r t 4 - 0 6 ) = C 1Q cos (w r t+ 0 4 )+ C lx cos(w r t+ 0 5 ) 


D 5 = 4M K x w r (a -HD r -w m )/(3E 1 ) 
D 6 « 8S a q u>_ U m /< 3E P 


(5.60) 
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If equation (5.57) is to give a constant value, wo should 
satisfy following relations 

l[D 5 cos(u> t+a) + sin(w ra t+a) ] = O (5.61) 

M[D 5 cos(2a> m t+2a) + sin(2a> m t+2ct) ] =0 (5,62) 

MC 10 = 0 (5.63) 

IC i2 = 0 (5.64) 

equation (5.62) can be- written as 

hVf ( D^+D^ ) cos(2w m t+2a+0 7 ) = 0 (5.65) 

v/hcre 

tan 0 7 = -(D^/Dc,) (5.67) 

Equation (5.65) is satisfied only if M = 0 or and are 
zero (5,68) 

M = 0 implies the case of no modulation and is not of interest. 
Thus , wo should have 

— 0 (5,70) 

D 6 =0 (5.71) 

and from equation (5.63) 

C 12 


= 0 


(5.72) 
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Equations (o.70) ato (5.72) arc the only possible condition 
which satisfy equation (5.61) to (5,64), for cosinusoidal 
modulation. Substituting these in equation (5,57) gives 

T q = 0 (5.73) 

This corresponds to the trivial; case of no stator current. 

Thus it is not possible to reduce current harmonics to zero, 
using the cosinusoidal modulation. 

5.4.2 Frequency domain analysis for the cosinusoidal 
modulation 

The effects on the torque spectrum due to the variations 
in the amplitude, frequency and the phase of the cosinusoidal 
modulating waveform have been studied in this section. The 
torque spectrum produced by the assumed stator current is 
computed using the method of Section 4.3.2. 

The torque spectrum for the cases with w , of equation 
(5,38), equal to 3,6,12 and 18 times the inverter frequency have 
boon studied. When s tudying the nature of torque spectrum it 
is found that for a particular choice w m and there is a value 
of a which gives the maximum performance index • Also 

there is a optimum combination of amplitude and a frequency 
which gives the lowest torque harmonics. 


The results have been summarised as follows : 
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(i) = 6w ; Table 5.2 gives the 6th and the average 
torque for various amplitudes and phases for the case of or, 
equal to 6to and a fixed constant rotor speed. This shows that 

(a) There is a distinct phase angle for a every choice of 
amplitude which gives the maximum performance index for 6th 
harmonic . 

(b) There is a optimum amplitude which for a particular phase 
angle gives zero sixth harmonic torque. 

(c) As the amplitude of modulating waveform increases the 
phase angle at which the minimum 6th harmonic torque is obtained 
decreases . 

(d) The variations in 'the average torque produced duo to 
modulation are not much significant. 

(ii) In Table 5.3, the effects on 12th harmonic torque for the 
case of a) = 6m has been taken up. It can be noted that we are 
able to obtain the reduction in the 12th harmonic torque with 
modulating current frequency equal to six times the inverter 
frequency. In this case also the property of the optimum phase 
and amplitude has been observed. 

(iii) co m = 3w i Table 5.4 lists the results for u) m equal to 
three times the inverter frequency. This typo of harmonic can 
be present in the invertor input dc current due to imperfect 
filtering. It is noted that there is an improvement In the 
ratio of average to 6th harmonic torque and the feature of 
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optimum amplitude and phase angle is present and is similar 
case (i). But for this case the change in the average torq\ 
is found to be significant. Even for a particular amplitude 
of the modulating waveform, there is a significant variation 
in the average torque with the variations in the phase angle 

(iv) = 12m ; The reduction in 6th harmonic torque by 

modulating waveform with equal to twelve times the invorti 
frequency has also been studied. Table 5.5 lists the result 
obtained for this analysis. It is observed that in this case 
the variation in the sixth harmonic torque due to the varia- 
tions in the amplitude .and the phase of the modulating wavefo 
is small, as compared to the frequency of 6w. Also, hero a 
1 arge scale modulation is required for better performance. He 
no optimum property with respect to the amplitude of modulatii 
current is observed in the possible range of amplitude of modi 
1 a tion, 

(v) Table 5.6, lists the result for the reduction of 6th har- 
monic torque with = 18o). Here again the variations in 6th 
harmonic torques are observed to be small. Those are oven 
smaller then for the case (iv) above. Here also no optimum 
property is observed. I 

Thus it can bo said that in case the cosinusoidal modula- 
tion is being used, the modulating waveform of the frequency 
six times the inverter frequency gives the best control for 
obtaining tho desired optimum performance. 
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Fig 5.3: Performance index Vs. rotor speed 
with p -- - a and variable K 
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H ARMfjR XF 
torque 

Tav/Thdr 

0.000 

14.64299 

4.08537 

3.58425 

0.000 

10.0 0 ( j 

72.000 

1 OR. 000 

180. 000 

216.000 
252.000 

i y |. ooo 

.’24. uOO 

14.53567 
14,55908 
14.61895 
14.68854 
14.74123 
14. 76073 
14.74200 
14.66987 
14.620J7 
14.56001 

4,30447 

3.07431 

2.21631 

2.52967 

3,57657 

4.6420O 

5.46017 

5.90231 

S. 87035 
5.31704 

3.37689 

4 . 7 j r 7 3 

6.59607 

5.80650 

4.12101 

3.17982 

2.69991 

2.4H883 

2 . 4 y 0 5 4 
7.73837 

0.000 

36.000 

72.000 
108.000 

144. 000 
16* .000 

216.000 
252. Off) 
28 8 . 000 
324.000 

14.43877 

14.49025 

14.61773 

14.75703 

14.85474 

14.88888 

14.85610 

14.75956 

14.62069 

.14.49229 

5.19065 

2.90330 

0.41514 

2.07946 

4.13491 

5.76392 

6.96377 

7.66192 

7 . 6 H 1 0 0 
6.85735 

2.78109 

4.9yQ8? 

36.21104 

7,09659 

3.59257 

2.58312 

2.13334 

1.92635 

1 .90 349 
2.11340 

J.000 

36.000 

72.000 
108,000 ■ 

144.000 

180.000 
216.000 

252.000 

288.000 
324, OOO 

14.35228 

14.43652 

14.63934 

14.84846 

14.98351 

15.02744 

14.90529 

14.85210 

14.64394 

14.43984 

6.46640 
3.73b84 
1.77840 
3.35017 
5.49003 
7.19967 
8,50947 
9,36304 
9.4661 8 
6.51 164 

2.21952 

3. 8o330 
8.23175 
4.43215 
7.72922 
2.08724 
1.76101 

1 .58625 
1.54697 
1,69648 

76.500 

76,600 

76.700 

76.800 

14,63818 
14.63861 
.14.6390 4 
14,63946 

0,00724 

0.00097 

0.00775 

0.01440 

.im:3 
2018 *0 
1016,4 



TABLE 5.3 

1 1 ***** **** 


FUR REPuC'f I ;JM r, J- 12 HAR, * IXH MOD, FREQ. AS 
6 TIMES MV. FREQ. A 1’ ROTHP SPEFD=52y . RPh 
n I'li XDC=b.O A VP 


AMP OF 

ANGLE 'W 

average 

tiARMONjC 

TavFTPar 

TOP (AMP ) 

MOD (DEG) 

TP ROUE 

TORQUE 

0.000 

0,000 

14.64299 

1,98144 

7.4 

1.000 

O.Ot, 0 

M. 4 3877 

2.o3fa43 

5.5 

1.000 

30.000 

14.49025 

2.48176 

5.9 

1.000 

72.000 

14.61773 

2.16417 

6.9. 

1.000 

108.000 

14.75703 

1.78471 

9.3 

1,000 

144.000 

14.85474 

1.45322 

10.2 

1 .000 

1 B 0,0 00 

14,88888 

i. 41 75 J 

It. 3 

t.000 

216.010} 

14,85610 

1.46165 

10.2 

1.000 

252.000 

14.75956 

1.83353 

B.O 

1.000 

2b 3. 000 

14.62069 

2.25776 

6.5 

1 .000 

324.000 

14.49229 

2.55982 

5.7 

2.000 

0.000 

1 4.27621 

3.29186 

4.4 

2.000 

36.000 

11.39788 

2.95252 

4,9 

2.000 

7 2 . OOO 

14,68378 

2 . 35866 

6.2 

2.000 

1 OB. 000 

14,96281 

1 ,67168 

9.0 

2.000 

144.000 

15.12765 

1.03165 

14.7 

2.000 

18.0,0 00 

15.17643 

0 . o 4 8 2 8 

23.4 

9.»0" 

216.000 

15.12957 

0,95773 

15.3 

2.000 

257.000 

14.96746 

1.75221 

8.5 

l:M 

283.000 

14.69012 

2.60522 

5.6 

324.000 

14,40265 • 

3,17808 

4.5 

3.000 

0.000 

14.15532 

3.91750 

3.6 

3.000 

36.000 

14.36588 

3,38989 

4.2 

3.000 

72.000 

1 1,84115 

2.54987 

5.8 

3.000 

1 OB . OOO 

15.26031 

l,o3751 

9.3 

3.000 

Ml. UNO 

15.46141 

0.80943 

19.1 

3,000 

1 S3 L. 000 

15.50565 

0.14185 

109.3 

3,000 

216,000 

15.46341 

0.63600 

24*3 

3,000 

252.000 

15,26667 

1.77166 

8.6 

3.000 

2 80, 000 

14.85130 

3.01507 

4.9 

3.000 

32 1.000 

14.37406 

3.82979 

3.8 

4,000 

0,000 

14.07609 

4.54330 

3.1 

4,000 

36,000 

14.39424 

3*79217 

3.9 

4,000 . 

4,000 

1 2.0 >0 
10 §,000 

15.08983 

15,64959 

2.72842 

1166758 

5.5 

9.4 

4.000 

4.000 

i •’ 4 . • i 

180,000 


0*90505 

0*75962 

17.5 

20.9 

4, 

4.000 

71 6.000 

252.000 

1 : 11 ? 

0,86196 

1*91391 

18.1 

8,2 

4,00m 

2d 3.0 00 

15,10421 

3.4P007 

4*3 

4.000 

324,000 

14,40652 

4.51192 

3.2 



1 { iJ ! * '* 6 FAR, 

3 J'l lS i , /. fpf.q. j\ r 
WITH IDC=6.0 


WITH 
HO TOP. 
A IPS 


MOD . FPEO, 
SPFEDs520, 


AS 

PPM 


A, if’ OF 
MOO (AMP) 

0.000 

0.500 
0.500 
0.500 
0.500 
0.500 
0.500 
0.500 
0.5 Op 
0.500 
0.500 

1.000 

1.000 

1.000 

1 . 000 

1.000 

1 ,000 

1.000 

1.000 
1*000 

1.000 

2.000 

2.000 

2.000 

2.000 

2.000 

2,000 

2.000 

2.000 

2.000 

2.000 


AUGUF .10 
MOD (DUG ) 


AVER ATE 
'frjpgJJE 


,0'iO 11.64299 


(.,000 

36.000 
72.00 < 

108.000 

lbo.opo 

210,000 

252.000 

288.000 

324.000 

0,0 0 0 

36.000 

72.000 

108.000 

144.000 

180.000 

216.000 

252.000 

288.000 

324.000 

. 0.000 
36.0(50 
72. OOP 

108.000 

1 44.000 

180.000 

216.000 

252.000 

288.000 

324.000 


14.67053 

13.29620 

12.47106 

12.4o395 

13.27705 

14.64560 

16.07634 

16.99415 

17,00232 

16.0972? 

14.72843 

12,03599 

10,47836 

10.46519 

11.99939 

14,67897 

17.59626 

19.52454 

19.54195 

17.63972 


14.93532 
9.77526 
“ 03066 
00855 
.70892 
14,83639" 
20.89581 
12301 
16206 


7 . 

7, 

9, 


25 

25 


20.98959 


H APHONIC 
TOP 0 OF 

4.08537 

5 . 66336 
4.93793 
3.946PJ 
2.99474 
2.42174 
2.55542 
3.33727 
4. <4 1 237 
5.3741 7 
5.84672 

7.21695 
5.74623 
3.92273 
2.15426 
0.85713 
1 .48910 
3,01649 
4.97940 
0.83709 
7.70583 


16.18541 
7.01585 
3.93735 
1.58291 
1.99219 
3,71072 
4.68565 
6.67766 
10.24406 
3 1.71798 


Tav/Tnar 


3.6 

2.6 

2.7 

3.2 

4.2 

K rr 

5 .*7 

4.6 

3.9 

3.2 

2.8 

2.0 

2.1 

2*7 

4.9 
14.0 

9.9 

5.8 

3.9 

2.9 

2.3 

1,5 

1.4 
1.8 

4.4 

4.9 
4.0 

4.5 

3.7 

2.5 

1.8 



.TABLE 5,5 

*********** 


for REDunrofi 

12 TIMES X w V , 

* 1 1 f H 


»F b HAK. 
FREQ. AT 

II)C=6AMPS 


WITH 

ROTOR 


AFP OF 
Ml *0 ( A IP ) 

0.000 

0,500 
0,500 
0.500 
U.,500 
0.500 
0. 500 
0.500 
0.500 

S # 5 0 0 
1.500 
.1.500 
1.500 

1 .50c 
1.500 
1.500 
1.500 


*500 
.500 
. 5 0 0 
,500 
, 500 
,500 
,500 
,500 


4.000 

4.000 
4.000 
4.000 
4.000 
4.000 

4.000 

4.000 

5.000 

5.000 

5.000 

5.000 

5.000 
5,00(i 

5.000 

6.000 
6.000 
6.000 
f.OOo 
6,000 
6,000 


A .COE It 
HOO(OSG) 


average 

TORyllE 


0.000 14,64209 


0.000 
36,000 
72 1 6 00 
1 OR. 000 
t yo. ooo 

252.000 

? w w . o o o 

324.000 

0.000 

36.000 

72.000 

1 OR. 000 

180.000 

252.000 

268.000 

324.000 

0.000 

36.000 
, 72,000 

100. 000 

180.000 

252.000 

288.000 

324.000 

P.000 
36.000 
. 7 2.000 

1 08.000 
180,000 

252.000 

288.000 

324.000 

0.000 

36.000 

72.000 

108.000 
180.000 

25 2.000 

324.000 


i 

1 

14 
14 
14 
1.4 
14 
1 4 


.61425 

.62053 

.63702 

.65629 

.67650 

.65697 

.63776 

.62106 

*57111 
.59369 
.6491 1 
.70707 
,75787 
,70898 

• 65152 

• 59542 


1 1.54709 
1 1.59072 
1 1.69325 
14,79014 
14.85835 
14.79303 
14.69754 
14.59406 

14.54690 

14.63111 

14,81956 

14.97523 

15.04492 

14.97920 

14.82707 

14.63753 

JH7066 
14.68793 
14.94383 
,13897 

urn 


IP' 


6.000 


14.69685 

0,0 to ! A. 61 364 

ifMui 

7 |, 0 C 0 15.10014 

^ *000 15.33493 

M:S<w ■-•' se5? 

124.000 


Mon. F»EO, 
SPEEns526, 


H ARMOR TC 
TORQUE 


AS 

RPM 


4. 18537 


3.92203 
3.99079 
4.0964(1 
19314 
25071 
10386 
99823 
92519 


3.60045 
3 . 65776 
4.19937 
4.45898 
4.53383 
4.23337 
3,69725 
3,61843 

3.28704 
3.79980 
4; 39372 
4.77413 
4,92205 
4.48632 
3.91029 
3.33770 

2.83720 

3.85196 

4.80953 

5.30255 

5.43660 

5.07838 

4.16614 

2.99266 

2.55568 

:?5?n 

|:fp77 

Uflii 


9462 


H-jfm 



Tav/r^ar 


3 . 6 

3.7 
3.7 
3.6 
3.5 

3.5 

3.6 

3.7 

3.7 

4.0 

3.8 

3.5 

3.3 

3.2 

3.5 

3,s 

4.0 


4 

3 

3 

3 

3 

3 

3 

4 

5 
3 
3 
2 
2 
2 

3 

4 


4 

3 

3 

1 

0 

3 
8 

4 

1 

8 

I 

8 

8 

9 

6 

6 


5.7 

3.7 
2.9 

2.7 
2.6 
2.7 
5.2 

6.4 

3.6 

2.7 

2.5 

2.5 

2.6 
5.4 


TABLE* 5.6 

*********** 


FUR 

18 


pe nuryr i on 
times T *! V 


OF 6 HAR. WITH MUD. FPEO 
. FREQ. A T ROTUR SPFEn=52M 
WITH [ DC-6 . 0 AmPS 


AS 

Rpr 


AMP OF 

ANGLE OF 

average 

HARMON TC 

Tav/Tnar 

MOD ( AMP) 

w.JDCnEG) 

TORQUE 

TO ROUE 


O.noo 

0.000 

14.64299 

4.08537 

3.6 

0.500 

0.000 

14,62872 

4.02507 

3.6 

• '.500 

3 6 , 0 u 0 

14.63185 

4.04398 

3.6 

0.500 

72.000 

14.64013 

4.07950 

3 .6 

0.500 

108.000 

14.64'»91 

4.11604 

3.6 

0.500 

144,000 

14.65743 

4.14048 

3.5 

0,500 

1 60.000 

1 4.66030 

4.1 4600 

3.5 

0.500 

216,000 

14.65771 

4.13130 

3.5 

0.500 

252.000 

14.65039 

4.09y9t> 

3.6 

0.500 

268.000 

14.64065 

4.06181 

3.6 

0.500 

324.000 

14.63219 

4.03218 

3.6 

1 .500 

0,000 

14.60927 

3.90550 

3.7 

2.500 

0.000 

14.60194 

3.78751 

3.9 

3.500 

0.000 

14.60675 

3.67132 

4.0 

5.000 

0.000 

14.63670 

3.50100 

4.2 

6.000 

0.000 

14.67183 

3.39059 

4.3 

6.000 

36.000 

14,74586 

3.9135U 

3,8 

6.000 

72.000 

14.90824 

4.51036 

3.3 

6.000 

100.01 0 

15.02815 

4.81282 

3.1 

6.000 

144.000 

15.05698 

4.84931 

3.1 

6.000 

180.000 

.16.050 79 

4,82455 

3.1 

6.000 

216.000 

15.05619 

4189362 

3.1 

6.000 

252.000 

15.031,39 

4.83267 

-3.1 

6.000 

6.000 

280.000 

324.000 

14.91709 

14.75113 

4.33743 

3 .5°97 4 

3.4 

4.1 


6,1 


CHAPTER 6 

CALCULATION OF CURRENT PROFILE FOR SPECIFIED 

TCRqUE WAVEFORM 


6.1 INTRODUCTION 

In the previous chapters the form of the stator current 
of the motor has been assumed and then the torque waveform 
has been obtained. In this chapter, the inverse problem of 
the computation of the stator current to obtain a assumed 
torque waveform has been studied. This stator current is 
assumc-G to bo produced by an ideal throe phase current source 
inverter whoso input dc current is being modulated, 

A nonlinear second order equation has been derived for an 
interval of the inverter. The solution of this equation gives 
the stator cunont wave f orm for the assumed interval. The 
stator current waveform needs to bo computed only for one 
interval as the invertor input current needs to be modulated 
at six times the inverter frequency (Sec. 5.2). Thus the 
scar or current during complete inverter cycle can be obtained 
from t\S‘ solution of stator current during any one interval. 

The equation to be solved Is of the second order. The 
solution of this equation requires two boundary conditions. 
Those have been obtained, using the fact that under steady 
spate, xhu rotor mmf must move by 60° in space in each T/6 
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interval so that each interval is similar. 

In this chapter, the stator current waveform required to 
obtain ? constant torque has been derived. For this case, 
the nonlinear second order equation reduces to two first 
order nonlinear equations. The analytical solution of 
these equations for the stationary rotor case has been 
obtained. This gives a exponential modulation as obtained 
in Section 5.3.1. For the rotating rotor case, those equa- 
tions have been solved through two techniques, namely phase 
plane analysis and numerical analysis. The numerical solu- 
tion has been done by a computer program. The solutions for 
the stator current which give the constant torque have been 
thus obtained for various values of the rotor speeds. 

In this chapter the interval of the invertor for which 
the analysis has been done is interval I (Fig. 2.1). It 
should again be noted here that the inverter has been assumed 
to bo ideal for the analysis. 

6.2 DERIVATION OF EQUATIONS FOR STATOR CUR IENT TO PRODUCE 

DESIRED TORQUE ' JAVEFORM 

In this section, th_ equation for stator current wave- 
form during interval I is obtained, which gives an arbitary 
torque waveform T (t). This equation is a second order non- 
linear equations. The two boundary value conditions required 
to solve this equation have also been obtained in this section* 
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6,2.1 Derivations of nonlinear equation 

Let us assume a General form i„(t) of the ’ a ! phase 

a 

current during interval I, which gives the desired torque 

waveform T (t). During interval I. for the ideal inverter 

q 


x b “ x a 


i c = 0 

(wo) 

Thus from equations the dq stator currents can be 

obtained as 


( 6 . 1 ) 

( 6 . 2 ) 


i dl = i a (t) 


V 


V3 


(6.3) 

(6.4) 
, 33 

Substituting those values of i^ and i ^ in equation (4.4^7? 
we obtain 

X, 


T ( t) = “i'Li.(t) [X-, + ™ ] 
q c a 1 f 3 


(6.5) 


In equation (6.5), T (t) is an assumed function of time and 

Q 

i (t) has to be evaluated, 
a 

From equations (2.52) and (2.53) 


( ? 2 +2ap+3 2 +.o2)x 1 = K x ai ql + K 1 u> r i dl +K 1 pi ql 


(p 2 +2a?+a 2 +f ) 2 )X 2 = *^ 1 a io 1 ~ K l t ’r i ql +K l pi 


dl 


( 6 . 6 ) 

(6.7) 
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where to^ is rotor speed and a, are given by equa 
ana (2.53a) respectively, 

I i om equa Lions (6.6) and (6.7), we have 
(p 2 +2 ap+a 2 +w 2)(x i + £§) = K ia (i ql + ^|1) + 


tion (2.18) 


K l p( V + K l“r ti dl - 


Substituting for i^ and i ^ in above equation from equations 
(6.3) and (6.4), we have 


(p 2 +2ap+a 2 +o?)( Xl + n|) = | K x w r i a 


( 6 . 8 ) 


Substituting for i (t) in above equation from equation (6.5) 

ct 

nives 


(p 2 +2ap+a 2 -Ho 2 )(X 1 + 4) 


4K i u r yiL 

. { - Xj 

( x l + TS } 


( 6 . 9 ) 


This can be written as 


( p 2 -f-2ap+a 2 +w 2 ) (y(t) ) 


( 4K 1% T a (t) 

'' "yW 


(6.10) 


where y( t) = (Xj^ + y^|) 


( 6 . 11 ) 


Equation (6.5) using equation (6.11) gives 
i a (t) = -T (t)/(M c y(t) ) 


( 6 . 12 ) 
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Equation (6,10) can be solved to obtain y(t). This on 
substitution in (6.12) gives the form of stator current 
which gives the desired torque waveform T (t)» Equation (6,10 
being of order two requires two boundary value conditions. 
These have been obtained in Sec. 6.2.2. 

6.2.2 Derivation of the boundary conditions 

The values of dq currents at the end of the inverter 
period are related to the initial values by the equation 
(2.77) as 


X 2 (T c ) * (1/2) X^o) - (V3/2) X 1 (o) (6.12) 

X X (T C ) = (f 3/2)X2 (o) + (1/2) X 1 (o) (6.13) 

where T_ is the time of interval I, that is, (u/3to), to being 
tho inverter frequency. Thus, Xj(T c ) and X 2 (T C ) correspond 
to the pseudo rotor currents at the end of interval I and 
X^(o) and X 2 (o) at the start of the interval. These relations 
are used to compute two boundary value conditions in terms of 
y(t) . 

From equations (2.50) and (2.51) 

X i (o) = to r X 2 (o) - a X x (o) + K x i ql (o) 

X 2 (o) = "^r x i^°) “ 2 X 2 ( 0 ) + i dl (o) 


(6.14) 

(6.15) 
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Here 


X,( o) 


(6.16] 


Xr>( o) — 


?roQ equations (6.14) and (6.15) we have 


(6.17) 


( X-, ( o) + 


X 0 (o) 


X,(o) 


X 2 (o). 


-yg~) ~ ^(X^o) - -p—)-a(X 1 (o) + ) + 


j l ( o ) 

+ K l^ql (o) + T3~" ) 

Using definition of y(t) from equation (6.11) and equations 
(6.3) and (6.4) in above, we obtain 


Y ( o) = o> r (X2(o) 


X-j ( o) 


-)-ay(o) 


(6.18) 


y(o) = 


(6.19) 


From equation (6.12) and (6.13) 

X 0 (T ) i X 0 ( o) lf . X, ( o) 

^ X l^ T c )+ ~T3 ) = 2 ^ x i(°) + ~T T“ > + *2 ( x 2 ( ° ) - “Y3~ ) 


Using the definition of y(t) from equation (6.11) and substi- 

X,(o) 

tuting for (X 2 (o) - “-y^— -) from equation (6.18) in above we 
obtain 


y(T ) = 1/2(1+ ^) y(o) + &- y(o) 


(6.20) 
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vvh 3 


ore 


Y( T ) = y( t=T ) 


( 6 . 21 ) 


similarly i arorn equations (2.50) and (2.51) 


G / 

Ut ^ X l + 


n(T c ) 


,Xi+p> | t=T =w r (X 1 (T c ) ~ + 


C 2CT c ) 


0 + K,(i ,(T r ) + 


1 dl^ T c^ 


Using equations (6.11), (6,3) and (6.4) in above, wc obtain 


X ( T ) 

( i c ) = o> r (>C^- > ( T r ) - -ij---- 0 --) - ayCTJ) 


(6.22) 


whore 


y(T c ) 


«Ly(t) j 

dt t=I\ 


(6.23) 


From equation -s (6.2) and (6.13) 


X,( TJ , X.(o) 

(5C 2 (T c ) -4s— £-) = 1/2(X2 (o) ^3 ) 


-if O Xn(o) 

'f (X x (0) + -yg ) 


Using oquatioms (6.11) and (6.18) this gives 

(X (T )- -) = — — ~ (y( o)+ay( o) ) - y(o) (6.24) 

2 c' M (2u>_) 

j. 


Substituting ir ncorn equation (6.24) into equation (6.22) we 


o btain 


Y(x ) = i (i+ y(o) - w r (i+ \) y(o) 

C r* ti) 


(6.25) 
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The equations (6.20) and (6.25) gives the desired boundar' 
value conditions required to solve equation (6.10). These 
relate the end point value and slope of variable y(t) to the 
initial value and slope of this variable. 

6.2.3 Derivation of the first order equations relating y and 
(dy/dt) 

Let a new variable z be defined as 


py = dy/dt = z 
Then this implies that 


P 2y = (dy/dt) = dz/dt =(dz/dy) (dy/dt) 

That is. 


P 


2 


y = 


Z 


dz 

dy 


(6.27) 


(6.28) 


Substituting equations (6.27) and (6.28) into equation (6.9) 
we got 


dz 

dy 


+ 2az 


4 

3 


K.o) T (t) , 0 0 

— ~) ~ - (a 2 +to 2 )y 


(6.29) 


Thus the second order equation (6,9) can be reduced to two 
first order equations given by equations (6.27) and (6.28). The 
boundary value conditions of equations (6.20) and (6.25) become 


y(T c ) = \ y(o) + H z(o) 

r “ r 

2 

z(T ) = 2 (1 + z(o) - ^5 U (1+ %£) y( o) 


(6.30) 

(6.31) 
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6.2.4 General properties of the solution of the derived 
equations 

There are certain properties of the solution y(t) which 
can be directly noted from the derived differential equations 
(6.10) and (6.29) 

(i) It can be seen from equation (6.10) that if y(t) is a 
solution for torque T Q (t), then (F) y(t) is also a 
solution for the torque value of F.T (t). Therefore, from 
equation (6.12) the current should be increased by a factor 
of (F) 1 ^ 2 to increase the torque by a factor F. 

(ii) If z = f(y) is a solution of equation (6.29) then from 
equation (6.29) it can bo seen that 

dz , _ dz 

<*Y |y 0 *z 0 dy |-Y 0 ,-z o 

This property eliminate the consideration of two quadrants 
of the phase plane plot [13] because of the symmetry* 

(iii) It is clear from equation (6.10) that if y(t) is a 
solution of this equation then -y(t) is also a solution, T^(t) 
is the desired torque waveform and has the frequency of six 
times the invertor frequency (See. 4. 2). Thus this property 
imply that if y(t) is a solution in interval I then ~y(t) is 
automatically the solution in interval IV. 
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(iv) Let us study the behaviour of equation (6.10) at the 
point whon the derivative of y(t) is zero, i.e., py = 0. At 
such a point the equation (6.10) reduces to 


( 2 2 2 , 

(p +a +u r )y 


4 K l“r T a (t) 
3 v(tT 


This implies that 


2 4 K l w r T q^ t ^ , 2 2s 

py = ’5 y(tj “ (a 


y(t) 


(6.32) 


During interval I, i (t) is always positive. 

Thus from equation (6.12), y(t) is negative during this interval. 
Thus equation (6.32) implies that d 2 y/dt 2 is always positive 
when dy/dt is zero. Thus, for assumed positive torque function 
t (t), if there is a possible solution of stator current, then 

* y M 

the form of juJ^) during interval I is one of the three types of 
those shown in Fig. 6.1. 

6.3 PHASE PLANE ANALYSIS FOR THE CASE OF CONSTANT TORQUE 

In this section the procedure to obtain the solution for 
tho stator current which produces a constant torque has been 
given. 

For the case of constant torque, T the equation (6.29) 
becomes 

dz _ „ A [SL + j 3 2 y ] -2 a 

dy z L y K 


(6.33) 
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where 


a 



K.u T 
l r gc 

w r 


(6.34a) 


P 2 = (a 2 +u 2 ) (6.34b) 

Fho equation (6.63) gives the relationship between y and its 
derivative z (= dy/dt). The phase plane plot [13] can, 
therefore, bo obtained from equation (6.33). 

The slope of the trajectories dz/dy = S say has following 
properties. 

(i) It tends to infinity as z tends to zero or y tends to 
zero. Wo note that from equation (6.34a) that a > o as we are 
dealing with the case of torque production due to rotor speeds 
loss then synchronous speed and hence T is positive. 

M. C 

(ii) If z‘is very largo cn -2a and is a constant independent 
of y. This is true therefore in all quadrants of y-z plane. 

(iii) Zero slope isoclin^fe. i.e. M = o 


§ = 0 « * = - h <7 + ^ < 6 - 35 > 

This defines a curve in the y-z plane. 

The locus of points in the phase plane wherein the trajectories 
have the same slope is referred to as isocline 
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This curve- has a maxima or minima when 


JL. (_ SL. + q~) = o 

2a ^ y 2 + ^ ; u 


This gives 


y = y 0 = + s" 


( 6.36) 


Corresponding to this y of equation (6.36) 


z = z. 


(2&LJL.) 1 


za y 


v = z = + VSL. o 
o — a P 


(6.37) 


(iv) For the isoclin with slope 3 , i.e 


d z/dy = S 


From equation (6.33), the equation of the curve in y~z plane for 


the isocline with slope is 


1 / o* * 0 2. v 

z = - — — (- + p y) 

(S+2a) y 


(6.38) 


i_ his curve has maxima or minima at 


c_ a. s 2. 


(3+2a) Y 


This gives 


(6.39) 
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This y m of equation (6.39) is same as y Q of equation (6.36). 
As equation (6.39) is independent of S, it implies that all the 
isocline curves have a minima at the same value of y given by- 
equation (6.39). From the equation (6.38) it is also clear 
that the isocline curves for different values of S is similar, 
except for a different constant multiple for z. 

(v) It can be seen from equation (6.33) that the reversing of 
the signs of z and y, we get the same equation. Thus, the 
phase plane plot need to be plotted only for 2 quadrants. 

(vi) From equation (6.33) 



For the zero slope isocline, i.e. for dz/dy =0, at the extreme 
point (y 0 ,z 0 ) from equation (6.36) y Q = + ^|r . Thus at this 
point (y z ) from equation (6.71) 

j2 

= 0 (6.41) 

dy^ 


Thus the point (y^z^) is a point of inflexion because at this 


point? 


dz 

dy 


and 


d^z 

— U are zero, 
dy 


Using the above information and S as a parameter it is possible 
to draw a sot of isoclines throughout the phase plane. These 
are then used to determine the trajectories graphically. The 
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method is shown in Fig. 6.2. The lines from point P, have 
slopes and S 2 which arc indicated by the short line segments 
crossing each isocline. They are extended until they meet the 
isocline indicating slope iv'^* The point P 2 is midway between 
the intersections. The process is then continued. The tra- 
jectory is plotted as a smooth curve through the points Pj,P 2 ,* 
which are thus determined. 


The trajectory satisfying the boundary value relations 

(6.30) and (6.31) gives us the desired stator current waveform. 
To find the stator current satisfying equations (6.30) and 

(6.31) following procedure is adopted. 

Let some arbitary value y(o) be assumed for solution. 

Then from equations (6.30) and (6.31), the equation for end 
point is given by 

z(t c ) * 2 (l + ^srHy(T c ) - *?(1+ ^§^)y(o)) + 

r r v 

2 

- ^ w r (l+ : o*r)y( o) 


X « G * $ 

z(T c ) = Y(T c ) + ivL, y(o) 


(6.42) 


where 


M 


1 





+ 



) 



^ ( i+ ^W) ] 


(6.43) 


(6.44) 
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Since y(o) is assumed constant, the equation (6.42) is a equa- 
tion of straight line. The intersection of this straight line 
with the trajectory passing through the point y(o) and z(o) 
gives one value of y(T c ). This is compared with the y(T ) 
obtained from equation (6,30). In case these are sane, the 
assumed z(o) for y(o) choice gives the solution else a new 
trajectory on the line y(o) = constant is followed. 

The phase plane plot with the trajectories for the case 
of Tq C = 4 Nt m and u r = 50 rpm has boon plotted in Fig. 6.3. 
Tnis has been plotted only for two quadrants because quadrants 
II and III are similar to quadrants IV and I respectively 
because of property (ii) (Sec. 6.2.4). 

To obtain the solution for different rotor speed, a 
different phase plane plot has to plotted. It can be shown 
that with the change of rotor speeds, wo can change the axis 
In such a manner that the isocline plots remain invariant. 
This can be seen as follows. 

Equation (6.38) gives the equation of isocline v/ith 
slope S. This can bo written as 

^ = “ Tsfe) ( y + (6.45) 



(6.46) 
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If the rotor speed is changed, the phase plane plot is now 

assumed for different torque T I given by 

qc new 


co 


’ qc new 


GO 


rf olrh 
■ r ^new 


(- 


o 

a +uo j 


r ‘new^ 




a +oo 


r 1 old 


'qc 1 old 


(6,47) 


with such a choice of Tq C , the right hand side of equation 
remains the same even for different rotor speed but with 
new value of constant torque. So, if the new z-axis scale is 
assumed as (z/a), the isocline plot will remain the same. It 
should be noted that since- the scale of z axis is changed, the 
value of S for each isocline will change but the plot of iso- 
clinnos need not be changed. 


The phase plane solution gives the values of y(t) as a 
function of y. The time interval can be evaluated from this 
information. The time interval, t-^ between y(o) and y^ is giver 
by 


*1 

t 1 = / (l/z)dy (6.48) 

y(o) 


Thus from this information of y(t) as a function of time, the 
stator current -waveform can bo evaluated. 


Since y(T ) point corresponds to the end of interval, 
thus from equation (6.48) 

y(t c ) 

T c = /. . (l/z)dy 
c y(o) 


(6.49) 
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A interval corresponds to a 60° interval. Thus inverter 
frequency, m, is given as 


w = it/3T 

c 


(6.50) 


Thus , 


TZ 


CO 


““yTTp™ - 

3[ / (l/z)dy] 

y(o) 


(6.51) 


6.4 GENERAL NUMERICAL SOLUTION FOR THE STATOR CURRENT 

In this section, the equation (6.10) has been solved to 
obtain the stator current waveform which produces the 
torque waveform T (t). For the case of stationary rotor, an 
analytical solution can be obtained for the stator current 
waveform which will produce the torque T^(t)» For the case 
of rotating rotor, the stator current has been obtained by the 
numerical solution of equations (6.29) and (6.27). These equa- 
tions are the first order equations corresponding to equation 
( 6 . 10 ) . 

6.4.1 Case of stationary rotor 

For the case of stationary rotor, (oj^ = O) , the equation 
(6.10) reduces as 

(p^+2ap+a^) y(t) =0 


(6.52) 
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For the solution of equation (6.52) these are two boundary 
value conditions. These have been obtained by substituting 
CD = C in equations (6.18) and (6.22) which gives 

y(o) = -ay(o) (6.53) 

y(T c )= -ay(T c ) (6.54) 

The solution of equation (6.52) is 

y(t) = Y 1 e wat + Y 2 t e~ at (6.55) 

where Y^ and are constants. These are evaluated using bound” 
ary value conditions. Substituting from equations (6.55) into 
(6.53) we obtain 

~aY x 4- Y 2 = -aY x 

i.o. f y 2 = 0 (6.56) 

r . . 

The equation (6.54) gives the condition simulatio n, as (6*56). 
Thus, from equation (6.12) 

MO = - ("HTd T q (t) e3t (6 - 57) 

cl 

The equation (6.57) gives the form of stator current which 
produces the torque waveform T (t). 

In case T (t) is assumed to be constant, it can be seen 
that equation (6.57) reduces to the form of stator phase current 
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obtained in Sgc. 5 ^ , 

* o.l.l, 

torque . 


which is shown to produce a constant 


6.4.2 Case of rotating rotor 

The solution of the equation (6.29) satisfying the 
boundary value coned itions given by equations (6.30) and (6.31) 
cano be obtained through numerical analysis. It should be 
nooed heie that hero T (t) is assumed to bo some- arbitary 

M 

desired function. 


Lot the variable y be y at some arbitary time t = t. 

0 o 

with z = z Q . Then from equation (6.29), 


dz 

dy 


y=y 0 



(« 2 ^)y 0 ] 


2a 


(6.58) 


To the first dogree of approximation, the values of time 
and variable z at y — (y Q +tiy), will be given as 




y=y G 


dy 


(6.59) 


t |y= (y c +Ay) - V Ay (6.60) 

where dy is a small change in y from the value of y Q . 

The stops for solving equation (6.29) satisfying equa- 
tions (6*30) and (6,31 ) f for a assumed constant rotor speed 
are given as follows. 
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» 

(i) Arbi t ary values of y(o) and z(c) are assumed for solving 
equation (6.39) 

(ii) y(T ), and z(l ), the end point values of the solution 

v v ' 

of equation (6.29) are calculated from equations (6.30) and 
(6.31) respectively, for the assumed values of y(o) and z(o). 

(iii) The equation (6.29) is solved from y(o) to y(T ) in stops 
of £ y, using equations (6.58) and (6.60) and the value of 

z(T ) is obtained. 

(iv) Values of z(T ) obtained in steps, (ii) and (iii) above 

V * 

are compared and the difference is calculated, 

(v) Stops (ii) to (v) are repeated with a modified value of 
z y keeping y Q as in step (i), till the difference in step 
(iv) is reduced below a limit. 

(vi) Step (iii) compute the time t, when the value of y has 
has become y(T c ). This time t corresponds to the value T c * 
From this, the inverter frequency, oo, can be computed, as 
from equation (6.50) 

w = tc/3T (6.61) 

The step (iii) gives the solution of equation (6.29) 

-when the assumed initial values of y(o) and z(o), give the 
difference of step (iv) less then a limit. From this, the 
stator current waveform during this interval can be obtained 
using equation (6.12). 
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Having obtained the stator current waveform for an 
interval of inverter, which is inverter I in our case, the 
stator current for other intervals can be obtained using the 
properties given in Sec. 5.2. 


This stator current will produce the assumed torque wave- 
form T ( t) . 

q 

6.5 NUMERICAL SOLUTION FOR THE CASE OF CONSTANT TORQUE 


Equation (6.29) has been solved for the case of constant 
torque. This equation for T^(t) = T CJC becomes 


„ dz , 0 „ 
2 Hy + 2az 


4 /^l w r T qc 
3 M_ 


) 


/ 2 2 \ 
(a -Hu )y 


(6.62) 


If a new variables are defined as 

Vr I * 

y * t * z = 

a 



(6.63) 


Then from equation (6.12), equation (6.62) becomes 


* dZ , o * 

z — ^ + 2az = 
dy 


4 , Wc , i. 
qc y 


, 2 ^ 2 , * 

- (a +a>_) y 


(6,64) 


* * 

The boundary conditions in terms of y and z are similar to 
equations (6,30) and (6.31) with y and z replaced by y and z 
respectively. 


This equation is solved using steps given in Section 6.4 
through a computer program, whose listing is given in the 
Appendix B at various values of rotor speeds* 
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Figs. 6.4 to Fig. 6.9 give the plots of initial points of 
solution and corresponding inverter frequency for various value 
of the rotor speeds. From these plots it can seen that 

(i) There is a lower limit of frequency in the plots. This is 
so because these plots have been drawn for a particular rotor 
speed and the constant torque is assumed to be positive. 

Thus the synchronous speed corresponding to the inverter fre- 
quency has always to be more then the rotor speed, 

(ii) The plot of inverter frequency with respect to the initial 

x 

value of the solution of y(t) is parabola. This parabola opens 
up as the rotor speed decreases. In case of parabolic form, 
there is a upper limit of frequency, 

(iii) There a lower limit on the value of y(o). This is so 
because protection of a constant value of torque has been 
assumed. The current is inversely proportional to y(t). The 
initial, slope of y(t) is negative and so the slope for the 
solution is always negative (Sec, 6.4.2 point (ii)). Thus the 
lower limit on y(o) corresponds such a initial current which 
will produce torque value greater then assumed torque value. 

Thus there is no solution during this region, 

(iv) For a particular value of inverter frequency there is only 
one solution for stator current which will give us the assumed 
constant torque. 
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(v) Points of interest of low slip one on the lower branch of 
the parabola. On this as the initial value y(o) value in the 
solution increases, z(o) value becomes more negative and the 
corresponding value of the inverter frequency decreases. 

(vi) From the Fig. 6,4 to Fig. 6.9 it is possible to plot the 
initial values (y(o), z(o)) for the solution for the various 
jof rotor speeds. This has been drawn for the inverter frequenc 
of 20 Hz in Fig. 6.10, for the constant torque at all rotor 
speeds. This plot can be used to obtain the stator current 
waveform to obtain constant torque at various rotor speeds for 
the inverter operating at certain constant frequency. 

Fig. 6.11 compares the solution of nonlinear equation 
at 500 rpm to produce a constant torque and the optimum exponen- 
tial waveform for this torque and rotor speed [sec. 5.3]. It 
is seen that there is a close correspondence. This implies that 
exponential waveform corresponds closely to the stator current 
required to produce constant torque. 
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speed = lOOrpm 
torque = 6 N-m 

y*= I (< a 




Fig 6-5; Initial Fbint Of The Solution 



\w / 





m * » * * * * i r\ / 



inverter freq(Hz) 


6 . 



y*1/i a 

speeds 500 r pm 
torque^ 30 N~m 





Fig kS • I nit iai Fbmt Of The Solution 



inverter freq (Hz) 


32 

28 



speed = 550rpm 
torques 30 N-m 

y=1/'a 



12 
4 " 


% +02 +04 +0-6 +0-8 +1-6 + 1-2 

no)— >• 




»n verier fre< 



300 ' 403 “" 500 ™ 600 
speed { r pm) - 




v «t i # p^/ 


6.33 


fqc = 22*1 

RPM = 500 
inverter freq^20 Hz 



n 4 - 


2 ■■ 

°0 " 1 5 ' 20 30 40 50 60 

wt ( in dogs }- — * 


Fig' 611 ; Current of 'a' phase to get constant torqui 
via nonlinear eq. £ min, sixtn harmonic 
•nr quo via exponent..-:,-' rr; 


** nN t i *«y 

, V.-% V h 



7.1 


CHAPTER 7 

CONCLUSIONS 

The analytical expressions for the referred rotor current; 
of the induction motor fed by ideal square wave current have 
been obtained. The rotor current has been computed by solving 
motor performance equations in dq frame. They are also comput- 
in';, by calculating the transfer function of the motor. It has 
been shown that both these methods give the same results. 

The analytical expressions for the induction motor being 
fed by stator current with cosinusoidal rise and fall during 
commutation have also been obtained. A more complex time 
domain expressions ere given. It has been shown that as 
to this solution approach the solution obtained for 

square wave current. 

The nature of the electromagnetic torque produced by 
the motor has boon studied. The general expression for the 
torque in a reference frame revolving at arbitary speed has 
been given. It has boon shown that the torque is independent 
of the speed of the frame. It has been shown that the torque 
harmonic have frequencies of the type 3pto, where u> is the 
frequency of stator current and p is a integer. For the case 
of motor being fed by an inverter, the torque frequencies are 
of the typo 6 poo. 
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The methods to compute electromagnetic torque have been 
given. This methods are applicable to both ideal and nonideal 
current source inverter. This torque has harmonics (Fig, 4,5). 
It has been shown that for computation of a torque harmonic of 
the order f 6p r , it is sufficient to consider the current har- 
monics of the order (6p-l) and (6p+l), (Table 4,2). Detailed 
calculation show that interaction of fundamental in rotor and 
those two harmonics in stator contribute the major percentage oj 
the torquo. 

The control of the torque harmonics by modulation of the dc 
input current has been studied in Chapters 5 and 6. It is showr 
that for the balanced system the modulating frequency should be 
six times the inverter frequency. The cases with exponential 
and cosinusoidal modulations have been studied to obtain con- 
stant torque. The study with exponential modulation shows % 

(i) It is possible to obtain a constant torque for stationary 
rotor with exponential modulation (Table 5.1). 

(ii) It is not possible to obtain constant torque for rotating 
rotor with exponential modulation, 

(iii) At speeds less then synchronous speed, the exponential 
modulation always improves the performance. At speeds greater 
then synchronous speed, there is deterioration in the perfor- 
mance with exponential modulation (Figs. 5.2 and 5.3). 
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(iv) At a given rotor speed the parameters of the exponentia 
modulation have a optimum value. With these values the sixth 
harmonic torque can be reduced to nearly z_zero (Fig. 5.4). 

The optimum choice varies with rctor speed (Fig. 5.5). 

"^ith cosinusoidal modulation, also, it is not possible to gel 

constant torque. Tables 5.2 to 5.6 show that 

(i) There is a distinct phase angle for every choice of 

amplitude and frequency of the cosinusoidal modulating wave- 
form which gives maximum performance index (T y/T^) 

(ii) = 6w, it is possible to reduce sixth harmonic torque 
nearly equal to zero (Table 5.2). 

(iii) There is not much variation in the average torque with 
a) = 6co, I2w or 18 ( 0 . With u> m = 3to there is a significant 
variation in the average torque. 

(iv) For u> m = 3(0 or 6w, the4e is an optimum amplitude for 
maximum performance index. For case of o> m = I2w or 18 ( 0 , 
this optimum nature is not observed. 

(v) There is not much improvement in performance index with 
ca = I2a> or 18w as compared to w n = 6w or 3w, 

(vi) It is possible to reduce 12th harmonic of the torque by 
current with oj = 6w (Table 5.3), 

iu 

(vii) It can be said that in case the cosinusoidal modulation 
is being used, the modulating waveform frequency o> m = 6w gives 


* 


best control. 
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The inverse problem cf computing the profile of the 
inverter output current t>: obtains the desired torque waveform 
has also been studied. The stator current for a desired torque 
waveform is given by the solution of a nonlinear second order 
equation. The expression for current to obtain a desired 
torque waveform at stationary rotor has been obtained. If the 
desired torque 'waveform is always positive, then the possible 
shapes of the stator current to obtain this torque is given 
by Fig. 6.1. The solution of the second order equation need 
two boundary value conditions. These have been obtained and 
these turn out to be independent of torque. This second order 
equation is reduced to two first order equations. These can 
bo solved by phase plane analysis or numerical integration. It 
can be- said that given the rotor speed and constant torque value 
it is possible to obtain a current waveshape which will give 
this performance. There is a range cf inverter frequency in 
'which this solution is possible (Fig. 6.4 to Fig. 6.9). If 
the invertor frequency is defined in this range then there is 
a unique current waveform. Outside this range, there is no 


solution. The solution of tho nonlinear equation gives current 
which is close to the exponential waveform (Fig. 6.11). 
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Details of tho slip ring induction 


Primary 
4C0 Volts 
5.0 Amperes 


3.0 H.P « 
50 Cycles 
1400 rpm 


The parameters of the motor are 

Stator leakage inductance L 11 

Rotor leakage inductance 

referred to the stator z 

Mutual inductance between 

stator and rotor phases 

referred to stator M 

Stator resistance 

Rotor resistance referred to 
stator 


motor used 

Secondary 
145 Volts 
10 Amperes 


= 0.2453 H 
= 0.2453 H 

= 0.2364 H 
= 1.6 Ohms 

= 3.31 Ohms 
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